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ABSTRACT 

We present mass models of a sample of 14 spiral and 14 SO galaxies that constrain their stellar 
and dark matter content. For each galaxy we derive the stellar mass distribution from near- 
infrared photometry under the assumptions of axisymmetry and a constant /T^-band stellar 
mass-to-light ratio, (M/L)/f^.. To this we add a dark halo assumed to follow a spherically sym- 
metric Navarro, Frenk & White (NEW) profile and a correlation between concentration and 
dark mass within the virial radius, Mdm- We solve the Jeans equations for the correspond- 
ing potential under the assumption of constant anisotropy in the meridional plane, j3,. By 
comparing the predicted second velocity moment to observed long-slit stellar kinematics, we 
determine the three best-fitting parameters of the model: (M/L)/f^, Mdm and jS,. These simple 
axisymmetric Jeans models are able to accurately reproduce the wide range of observed stellar 
kinematics, which typically extend to w 2-3 Re or, equivalently, w 0.5-1 /?25- We find a me- 
dian stellar mass-to-light ratio at A'^-band of I.Q9{M/L)ks.q with an rms scatter of 0.31. We 
present preliminary comparisons between this large sample of dynamically determined stellar 
mass-to-light ratios and the predictions of stellar population models. The stellar population 
models predict slightly lower mass-to-light ratios than we measure. The mass models contain 
a median of 15 per cent dark matter by mass within an effective radius R^. (defined here as the 
semi-major axis of the ellipse containing half the Ks-hand light), and 49 per cent within the 
optical radius R25. Dark and stellar matter contribute equally to the mass within a sphere of 
radius 4.1 R^. or 1.0/?25- There is no evidence of any significant difference in the dark matter 
content of the spirals and SOs in our sample. Although our sample contains barred galaxies, we 
argue a posteriori that the assumption of axisymmetry does not significantly affect our results. 
Models without dark matter are also able to satisfactorily reproduce the observed kinemat- 
ics in most cases. The improvement when a halo is added is statistically significant, however, 
and the stellar mass-to-light ratios of mass models with dark haloes match the independent 
expectations of stellar population models better. 

Key words: dark matter — galaxies: elliptical and lenticular, cD — galaxies: kinemat- 
ics and dynamics — galaxies: spiral — galaxies: stellar content — galaxies: fundamental pa- 
rameters 



1 INTRODUCTION 

In the dominant paradigm describing structure formation, initial 
fluctuations in density are enhanced by gravity until galaxies form 
in potential wells ( [White & Ree s 1978, Bl umenthal et al.||1984^ . 
Numerical simulations and semi-analytic models of galaxy forma- 
tion are able to reproduce many of the statistical characteristics of 
galaxies and some of their detailed features (e.g. Kauffmann et al. 
1993','Cole et al.|1994||Benson et al.|2003||Baugh|2006[|Somerville 
et al. 2008) . A detailed and consistent theory remains elusive, how- 
ever, and there are apparent contradictions between the predictions 
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of models and the observations (see, e.g. |Baugh|2006[ [Mayer et aL] 
|2008[ and references therein). 

The main stumbling blocks are the significant uncertainties 
and computational difficulties involved in capturing the baryonic 
physics that is crucial on sub-Mpc scales. Progress can be made in 
two ways. Firstly, using observational constraints which are thought 
to be insensitive to baryonic physics, one can circumvent these un- 
certainties and directly test the predictions of large-scle structure 
formation, which is in itself crucially important. Alternatively, one 
can test formation models with galaxy-scale observations. The out- 
comes of these tests can be used to refine and improve the models. 
Perhaps one of the most useful observational constraints for either 
approach is an understanding of the relative distribution of dark and 
luminous matter. In this work we aim to measure the radial dark 
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matter distribution in a sample of 28 edge-on disc galaxies. We 
make cosmologically-motivated assumptions, which allow us to lift 
certain modelling degeneracies. 

It is well-established that, within the optical disc, the kinemat- 
ics of high surface brightness spiral galaxies can be reproduced by 
maximal disc models, in which luminous material contributes the 
maximum amount consistent with the observed rotation curve (e.g. 
van Albada & Sancisi 1986, [Persic et al . 1996; Palunas & Wi Uiamsl 
20001. The stellar components of maximal disc models typically 
contribute 75-95 per cent of the rotational velocity at LlR^n^^, 
where i?disc is the scale length of the exponential disc and 2.2R^[gf. 
is the radius where the rotational velocity of the exponential disc 
peaks (,Sackett|1997^ . This permits dark haloes that comprise 10- 
45 per cent of t he total mass within 2.2 ^disc- In a Freeman disc 
I Freeman]|l970 R25, the radius of the 25 mag arcsec^^ isophote 
occurs at 3 iJdisc ■ 

Side-stepping the maximal disc assumption, Ratn am & Salucci| 
P000| l demonstrate that rotation curves can be adequately fitted 
without dark matter. [Kassin et aLlj2006^ avoided the maximal disc 
assumption entirely by using independent estimates of the stellar 
mass-to-light ratio inferred from a relationship with colour f Bell &| 
|de Jon g 2001 , B ell et al.|2003[ l. Most of their models were consis- 
tent with a maximal disc. Further evidence is provided by model- 
independent analysis of rotation curves at intermediate radii ( |Mc-| 
[Gaugh et al.|2007) . 

In the case of elliptical and SO galaxies, evidence from dynam- 
ical modelling and gravitational lensing studies suggests that dark 
matter makes only a small contribution to the total mass within i?e, 
the radius of the elliptical isophote containing half the light. For ex- 
ample, from the dynamical studies, [Gerhard et al.| ( |200T) find that 
10-40 per cent of the mass within Re is dark in their sample of 21 
ellipticals, |Borriello et al.|p003| l find 30 per cent dark matter within 
by constraining a sample of 221 ellipticals to lie on a fundamen- 
tal plane, Cappell ari et al.| ( |2006| l find a median of 30 per cent dark 
matter within R^ in their sample of 25 ellipticals and SOs, Thomas 
|et al.| ( (2007| find 10-50 per cent dark matter within Re for a sample 
of 17 ellipticals and SOs, and |Weijmans et al.|p008t find 55 per cent 
dark matter within 5 R^ in NGC 2974. Lensing studies find similar 
results. For example, |Rusin et al.| ( |2003| find 22 per cen t dark mat- 
ter within 2 R^ of 22 elliptical lenses, Koopmans et al.|(2006[ l find 
25 per cent dark matter within the Einstein radius of 15 elliptical 
lenses (the Einstein radius is approximately equal to the effective 
radius for galaxy scale lenses), and |Bolton et al.|(2008[ l find 38 per 
cent dark matter inside Rg for 53 elliptical lenses. 

However, the total absence of dark matter in ellipticals is 
sometimes excluded with only low significance due to intrinsic de- 
generacies in the dynamical models (see, e.g. |Romanowsky et ah] 
[2003] and the response by Dekel et al.|2005| i and the lack of kine- 
matic tracers at large radii. Moreover, in all cases R^ <C /?25- The 
constraints on both the mass of the halo and its extent are therefore 
less strong for ellipticals than for spirals. 

Perhaps the most compelling evidence for the dominance of 
stellar matter in either ellipticals or high surface brightness disc 
galaxies comes from the analysis of the bars and spiral arms often 
present in disc galaxies. These non-axisymmetric features lift the 
degeneracy between the contributions from luminous and dark mat- 
ter. This is done by assuming that dark matter is axisymmetric, so 
all non-circular motions can be attributed to the non-axisymmetric 
luminous component. This approach has been applied to both bars 
(e.g. jEnglmaier & Gerhard|[l999| [Weiner et af][200T| > and spiral 
arms (e.g. lKranz et al.^OOS ^ and provides results consistent with 
maximal disc studies. Further constraints come from W-body simu- 



lations of bars, which imply that a significant central dark compo- 
nent would slow or even destroy bars (e.g. jPebattista & Sellwoodj 
|2000^, while observation s systematically indicate that bars are fast 
(e.g. jAguerri et al.|2003[|Gerssen et al.|2003| l. 

To summarize, observational evidence indicates that the kine- 
matics of both spiral and elliptical galaxies can often be reproduced 
by mass models that include a sub-dominant contribution from dark 
matter in the optical or central regions. The specific amount re- 
quired and correlations between halo and luminous galaxy prop- 
erties are, however, unclear. Note that we have not mentioned 
low surface brightness and dwarf galaxies, which are dark matter- 
dominated at all radii (e.g. Persic et al.' 1996[|de Blok & McGaughj 
|1997[ |Verheijen|[l997| |SwaFers_ 1999j and may provide the most 
stringent constraints on halo shapes, once observational uncertain- 
ties are resolved (e.g. jde Blok et al.|2001||Swaters et al.|2003| and 
references therein). Such galaxies are, however, unrepresentative 
and extrapolating their results to systems with greater stellar masses 
may introduce biases. It is therefore crucial to place accurate con- 
straints on the dark matter content of giant, high surface brightness 
spiral, SO and elliptical galaxies. 

This work uses a modelling technique which is different but 
qualitatively similar to traditional mass decomposition and rotation 
curve analyses of spiral galaxies. The key differences are (i) the stel- 
lar components of our mass models are based on deep near-infrared 
photometry, which accurately traces the smooth stellar potential of 
the galaxy (see Section [2T^ , (ii) we lift a degeneracy by making as- 
sumptions about the shape of the dark halo that are motivated by the 
results of cosmological simulations and observational constraints 
(see Section [Z2t (iii) we lift a further degeneracy (and account for 
pressure support) by comparing the predicted second velocity mo- 
ment rather than the rotational velocity to the observed kinematics 
(see Section[23j. 

We apply this technique to a sample of 28 edge-on spiral and 
SO galaxies. We use edge-on galaxies because, under the assump- 
tion of axisymmetry, the deprojection is unique. Moreover, when a 
galaxy is close to edge-on (; = 90°), an inclination error does not 
propagate on to significant uncertainties in the kinematics in the 
plane of the galaxy (which is proportional to sini), and thus on to 
significant mass uncertainties. The mass models consist of a dark 
halo component and an unusually detailed parametrization of the 
projected light. Under justifiable assumptions, the mass model has 
three free parameters, the stellar mass-to-light ratio, the mass of 
the dark halo and the velocity anisotropy. We place constraints on 
these parameters by adjusting them so that the mass model predicts 
stellar kinematics that closely match those observed. We hope that 
these simple but powerful quantitative statements will be of use in 
constraining models of galaxy formation and evolution. 

In addition to the halo masses inferred, the constraints we 
place on the near-infrared stellar mass-to-light ratios are them- 
selves of great interest. There are significant difficulties in mod- 
elling the spectral energy distribution of stellar populations in the 
near-infrared, due to the importance of the complex thermally puls- 
ing asymptotic giant branch at these wavelengths (Maraston 2005 1. 
Moreover, the normalization of the models is uncertain due to a lack 
of knowledge about the initial mass function of stars. Dynamical 
measures of the near-infrared stellar mass-to-light ratio like ours, 
that do not depend on population models and seek to correctly ac- 
count for dark matter, therefore provide important independent tests 
of these models. 

A final goal of the study is to provide detailed constraints on 
dark matter in SO galaxies. The dark matter content of SOs is an 
important quantity to constrain because the dominant model of SO 
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formation as faded spirals predicts a simple and verifiable relation 
between the TuUy-Fisher relations of spirals and SOs (e.g. |Bedre-| 
|gal et al.|2006| and references therein). Galaxies that lie on a single 
Tully-Fisher relation linking their luminosities and rotational ve- 
locities are believed to have equal dynamical mass-to-light ratios. 
This will be investigated in more detail in a future paper (Williams, 
Bureau & Cappellari, in preperation). 

This paper is organised as follows: in Section |2] we describe 
the multi-Gaussian expansion used to model the luminous and dark 
mass distribution of the galaxies and the dark halo model adopted. 
We then give an overview of the Jeans modelling technique used to 
model the stellar kinematics. In Section[3]we present the sample of 
edge-on galaxies under consideration and describe the photometric 
and kinematic data. In Section|4]we present the results of our mod- 
elling of both the mass distribution and stellar kinematics of the 
sample, which we discuss in depth in Section [5] Finally, we sum- 
marise our conclusions and discuss possibilities for future work in 
Section|6] 



2 METHODS 

2.1 Modelling the luminous mass distribution 

We use multi-Gaussian expansions (MGEs) to create mass models 
that have a simple analytic form whose kinematics can easily be 
solved using the Jeans equations. MGE is a method of parametriz- 
ing an image of a galaxy as the sum of a finite number of two- 
dimensional Gaussian functions. The original application of two- 
dimensional Gaussians to galaxy images ( [Bendinelli|l"99T) was ex- 
tended to the non-circular case and to more general point spread 
functions (PSFs) by [Monnet et al.| ( [I992| > and [Emsellem et al.| 
( [1994[ l. Here we briefly summarize the method using a formalism 
due to |Cappellari| ( |2002l l, which we simplify for the special case 
of axisymmetric edge-on disc galaxies with a luminous component 
with a constant mass-to-light ratio. 

In a coordinate system (x, y, z), where x and z are centred 
on the galaxy nucleus and point along the major and minor axes 
in the plane of the sky while the y-axis points away from the ob- 
server, ILx z) the apparent surface brightness of a galaxy at an ar- 
bitrary waveband X can be written as a sum of A' two-dimensional 
Gaussians of apparent width CT; in the x-direction and d,-^, in the 
z-direction: 
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where L,- is the total luminosity of the ith Gaussian component. We 
further model the PSF as a circular Gaussian of width Opsf such that 
the intrinsic projected light distribution, deconvolved from seeing 
effects, is 
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where (7,-, the intrinsic width of the ;th Gaussian in the x-direction, 
and djqj, the intrinsic width in the z-direction, are given by 
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Once the projected light distribution is expressed in this sim- 
ple analytic form, it can be deprojected straightforwardly to give 



a full three-dimensional model of the light. In this work we de- 
project by assuming that the galaxy is axisymmetric, but other as- 
sumptions about the geometry are possible. Assuming axisymme- 
try and an edge-on view allows us to trivially transform from the 
(x, y, z) Cartesian coordinates of the Gaussians on the sky to the 
(R, (p, z) cylindrical system of the galaxy, where R is the galacto- 
centric radius, (j) the azimuthal angle and z the axis of symmetry 
of the galaxy. We also assume a constant stellar mass-to-light ra- 
tio (M/L)x to transform the stellar light at waveband A" to a stellar 
mass distribution. 

Together these assumptions of axisymmetry, an edge-on view 
and a constant mass-to-light ratio imply that the intrinsic mass dis- 
tribution of the luminous component of the galaxy can be written 
as 
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where a,- = L,/(v2;rcT/ ) ;/,■. 

In this work we determined an optimal MGE parametrization 
of the images described in Section|3]using the public fitting routine 
of|( 



Cappellari j2002|n which minimizes the quantity 
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for a given set of CT,-, qj and L,-, where M is the number of photo- 
metric data points Cj and £ is the apparent surface brightness of the 
model. 

A MGE is simply a sum of Gaussians reproducing the ob- 
served surface brightness. As such, the amplitude of each individual 
Gaussian, L/, is not constrained to be positive as long as the total 
luminosity and density are positive. Our sample consists of edge- 
on disc galaxies, which are particularly difficult to fit with solely 
positive Gaussian terms. As was demonstrated by [Bureau et al.| 
(20061, the major axis surface brightness profiles exhibit significant 
plateaus and secondary maxima. These cannot be fitted by a sum of 
concentric Gaussians with positive amplitudes, which is necessar- 
ily monotonically decreasing with radius. The complex rectilinear 
or concave two-dimensional structures visible in the isophotes of 
boxy and peanut-shaped (B/PS) bulges are similarly challenging. 

We therefore lifted the positivity constraint on L,- when mod- 
elling the luminous mass. In doing so, we encountered numerical 
issues with both the accuracy and stability of the fitting algorithm 
when using large numbers of Gaussians (~ 30). We solved these by 
reducing the maximum number of Gaussians in the sum to a rel- 
atively small number (~ 10), enabling double precision arithmetic 
to avoid cancellation errors, and finally, where necessary, tweaking 
the minimum surface brightness level down to which the fit was 
constrained by the photometry. We show a typical example of the 
improvement that is possible when terms with negative amplitude 
are allowed in Fig.[T] These improved mass models resulted in small 
but systematic improvements in the accuracy of the modelled kine- 
matics. That a more accurate description of the light gives a more 
accurate model of the kinematics gives us some confidence that our 
kinematic modelling methods and assumptions (described in Sec- 
tion |2.3[ l are not significantly flawed. 

The pixel-by-pixel absolute deviation from the photometry of 
our MGEs is typically 2-4 per cent. For a given constant {M/L)x 
these small errors are propogated linearly into the mass model. In 
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Figure 1. A comparison of the best MGE model with only positive am- 
plitude Gaussians (left-hand panels) and both positive and negative ampli- 
tude Gaussians (right-hand panels) for the same object, NGC 3957. The top 
plots show the A'5-band surface brightness of the image (filled contours) and 
model isophotes (solid lines). Contours are separated by 0.5 mag arcsec"^. 
The middle plots show the major axis surface brightness profiles of the im- 
age (points) and models (solid lines). The bottom panels show I'nns (points), 
the observed root mean square velocity, and ^12, the best-fitting second mo- 
ment found by solving the Jeans equations for the MGE model (solid lines). 
Note that allowing terms in the Gaussian expansion to have a negative am- 
plitude improves both the fit to the photometry and the accuracy with which 
the kinematics can be reproduced. In order to make the effect clear, the mass 
model shown here does not include a dark halo. 



most previous work, the observed surface brightness distribution 
of disc galaxies is parametrized using fits to azimuthally averaged 
radial profiles (e.g. [Ratnam & S alucci 2000[ [Kassin et al.||2006} 
or two-dimensional Sersic and/or exponential decompositions (e.g. 
[Gentile et "aLl|2004^ . While the terms in our MGEs lack any di- 
rect physical association with intrinsic components of the galaxies 
(bulge, disc, etc.), they reproduce the observed surface photometry 
more accurately than the simpler parametrizations, and their form 
is mathematically convenient for Jeans modelling (see Section|23J. 

At large radii, where the model is only weakly constrained by 
photometric data, models including negative Gaussians can look a 
little unphysical to the eye (see, e.g., NGC 1886 in Fig.[5](. Because 
there is so little light at these radii, however, we are confident that 
this does not affect the results significantly. We confirmed this by 
summing the light of all Gaussians in the MGE models, which can 
of course be computed analytically, to derive estimates of the total 
apparent magnitudes of the galaxies. As is shown in Fig. |2] these 
match the total apparent magnitudes presented in the 2MASS Ex- 
tended Source Catalog ( [Jarrett et al.|2000) to within 0.2 mag. This 
demonstrates that the excess light sometimes present in the outer 
regions of MGEs including negative terms is not significant, and 
that the photometric calibration described in Appendix [B] is reli- 
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Figure 2. Comparison of the apparent magnitudes of the sample galaxies 
computed by summing the light of all Gaussians in the best-fitting MGE 
models with the apparent magnitudes published in the 2MASS Extended 
Source Catalog. The 2MASS values are derived in the usual way using a 
growth curve. The methods typically agree to within 0.2 mag. SOs are shown 
as red circles and spirals as blue triangles. 

able. Fig.|2]also hints that the application of a growth curve method 
to the relatively shallow 2MASS images of faint objects results in 
total magnitudes that are systematically too faint (see also |Noorder-| 
|meer & Verheijen|2007) . 

2.2 Modelling the dark halo mass distribution 

If mass models are derived under the assumption that mass follows 
light then that implicitly assumes that dark matter is not a significant 
component by mass within the radius probed by the model. Strictly 
speaking, such models also admit the possibility that the spatial dis- 
tribution of dark matter closely matches that of the luminous matter, 
but we know of no physically motivated model of galaxy formation 
that predicts this behaviour. 

We extend our models by explicitly including a dark halo. Al- 
though the long-slit kinematic data for each galaxy typically reach 
~ 2~3Re or, equivalently, ~ 0.5-1 i?25! they are insufficient to al- 
low a us to constrain the shape of the halo of giant, high surface 
brightness galaxies like these (see, e.g. [Banerjee & Jog 2008). We 
therefore emphasise that we are not able to constrain the shapes of 
the dark haloes. Rather we attempt to probe the more limited ques- 
tion of what fraction by mass of each galaxy is dark assuming the 
haloes follow the particular one-parameter density profile described 
below. 

Our halo model is motivated by the results of A'-body simu- 
lations of cold dark matter (CDM) in an up-to-date cosmology. In 
such simulations the dark matter halo density Pdm is given by the 
spherically averaged NEW profile: 

P.S 



PDM(r) 



(r/r,)(I + ./.,)2 



(7) 



(e.g. |Navarroetal. 1997 , Bullock et al.'2001 :'Wechsler et al.'2002 1, 
were is a scale radius characterising the location of the 'break' in 
the profile and Ps is a corresponding inner density. 

Rewriting equation ^ in terms of Mdm, the total dark matter 
mass enclosed within the virial radius rvir, and defining a concen- 
tration parameter cvir = ''vir/rs implies 
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where Mdm = ^^PsrlA(c^\^) and 
A(x) = ln(l+x)- 



l 



(9) 



The virial radius is defined as the radius within which the mean 
density is Ap^rit and Pcrit.O = 

now. Throughout this work 
we adopt the cosmological parameters found by the Wilkinson Mi- 
crowave Anisotropy Probe five-year results (WMAP5), i.e. A = 95 . 1 
and the Hubble constant tfo = 100/zkms"' Mpc"' where /i = 0.70 
( |Komatsu et al.|2009| l. 

Following Napolitano et al. ( 2005), we use a key result from 
simulations of ACDM to eliminate one of the free parameters from 
equation l[8]l. The simulations demonstrate that concentration and 
halo mass are correlated: low mass haloes are more concentrated 
and the concentration-mass relation is well described by a single 
power law with a slope ~ —0.1 (e.g. [N avarro et al. 1997; Bullock 
let aL|200HrEke et al.|200H|Kuhlen et al..2005,,Neto et al..2007| l. 



We choose the fitting formula presented by Maccio et al. ( |2008[ l 
that is consistent with their WMAP5 simulations over the range 
1O'°M0 <Mvir < lO'^Ms: 



/i-11012M;^ 



-0.094 



(10) 



where is the total mass enclosed within the virial radius. 

We then rewrite this relationship in terms of the dark mass in- 
side the virial radius, rather than the total mass using the equation 
Mvir = a Mdm- A number of plausible choices are available to us for 
the constant a: to adopt the cosmological value, to neglect the con- 
tribution of baryons entirely, or to adopt some intermediate value. 
The cosmological value is derived from the ratios of matter density 
to critical density ilm = 0.258 and baryon density to critical density 
ilb = 0.0438, which together imply 

The most extreme 'missing baryon' scenario, which neglects the 
contribution of baryons to the total viral mass, implies a = 1. An in- 
termediate possibility is motivated by the results of observationally 
constrained halo occupation distribution methods and semi-analytic 
models, which imply that the stellar mass of a galaxy is typically 
around 3 per cent of that of its halo, i.e. a = 1.03. In this work we 
follow [Napolitano et al.|j2(305^ and adopt the cosmological value, 
a = 1 .20, which gives 



r(MDM) -9.195 



-0.094 



(12) 



However, we note that the choice of a makes almost no dif- 
ferent to our results. This is because a is immediately raised to the 
power —0.094 in order to define the concentration of the halo; a 
20 per cent change in a changes Cyjj at the 2 per cent level. If our 
choice is wrong and results in the introduction of such a small sys- 
tematic error in Cvir^ then the consequences for the parameters of 
the best-fitting mass models are in any case negligible compared to 
the observational errors. 

The dark halo density profile may therefore be written as 
a function of a single parameter, Mdm, using equations |8](, ^ 
and l |12^ . We perform a multi-Gaussian expansion of this one- 
dimensional, single parameter profile to allow us to easily include 
it in the potential for which we derive model kinematics. 

We refrain from including prescriptions for the effects of bary- 
onic contraction on our halo model (e.g. jBlumenthal et al.[|1986[ 
IGnedin et aL][2004l |Abadi et al.|[2009l >. We do this for two rea- 
sons. Firstly, our goal is not to determine the detailed shapes of dark 



haloes but rather the total dark and stellar masses in the optical parts 
of galaxies. Secondly, while our kinematic data do not allow us to 
constrain the halo shape, other observational evidence suggests that 
contracted NFW haloes do not reproduce observed kinematics, (e.g. 
[Gentile et al.|2004[[Kas"sin et al.|2006l [Thomas et al.|2007> . 

2.3 Modelling the stellar kinematics 

The most general dynamical methods are particle-based (e. g. [de 



'Lorenzi et al.'2Q07 ^ or orbit-based (e.g. ISchwarzschildl l979j Cap- 
pellari, et al. 2007| ]van den Bosch et aT 2008 ; Thomas e t al |2009^ . 
These powerful methods are so general, however, that a wide range 
of unrealistic models can be made to fit long-slit data, providing 
only weak constraints on the parameters of the mass models (see, 
e.g., fig. 2 of [Cappellari & McDermid|[2005l . In fact, the stellar 
kinematics of early-type fast-rotator galaxies are well-described by 
Jeans models with a cylindrically-aligned velocity ellipsoid with 
a constant flattening in the z-direction ( [Cappellari, et al.][2007[ 
[Thomas et al.|200^ . This has lead to the use of a Jeans modelling 
technique in which such a velocity ellipsoid is assumed. By varying 
the flattening of the velocity ellipsoid, the Jeans modelling approach 
has been shown to reproduce a wide range of two-dimensional ob- 
served kinematics in the central regions of early-type fast-rotators 
in the SAURON survey ( |Cappellari|2008l|Scott et al.|2009^ and ou t 
to 5i?e in the case of the edge-on SO NGC 2549 (| Weijmans|2009t . 
The large size of our sample (28 galaxies) allows us to test whether 
the assumptions provide a good description of our galaxies, but we 
note that our galaxies are at least as rotationally supported as the 
SAURON galaxies, so we expect to be even less vulnerable to as- 
sumptions about anisotropy. 

For convenience we provide an overview of the derivation and 
solution of the Jeans equations. This is essentially a summary of 
Sections 2 and 3.1 of'Cappel lari[j2008t . Under the assumptions re- 
quired for the coUisionless Boltzmann equation to hold (a smooth 
potential and steady state), and the further assumption of axisym- 
metry, the Jeans equations may be written in cylindrical coordinates 
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( |Jeans|I92"2j[Binney & Trem aine 2008 1. Here p is the density, * is 
the gravitational potential and we use the usual notation 



(15) 



Even if p and <I> are known (as is the case for our mass mod- 
els), equations \13\ and ( [14[ ) are still two equations with four un- 
knowns, v|, Vj, and vrv,, so they do not specify a unique so- 
lution. In order to close the equations one can assume a particular 
anisotropy, i.e. a relationship between the lengths of the axes of 
the velocity ellipsoid. We assume here that the velocity ellipsoid is 
aligned with the cylindrical coordinate system of the galaxy and that 
the anisotropy in the meridional plane is constant (v| = bv^ where 
is a constant). Under these assumptions, equations \13\ and l |14[ l 
reduce to: 



bP^'l-P^-l d{bpvl) 



dR 

dipvl) 



= -P 



dz 



d^ 

dR' 
d<i> 



(16) 
(17) 
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These equations are solved for the case of a density and potential 
described by a sum of Gaussians in I Cappellarij ( |2008 T|. Once the 
observable quantities in the solutions have been projected along the 
line-of-sight, a single integration gives a model of the second veloc- 
ity moment, ^2 (see equation [28] of |Cappellari|2008^ . We use the 
Jeans Anisotropic MGE (jam) routines to perform the calculation.' 
The second velocity moment is compared to the observed root mean 
square velocity, Vnns = (i'^ + C7^)'/^, where v is the observed line- 
of-sight velocity and a the observed line-of-sight velocity disper- 
sion. Vrnis is a good approximation to the true second moment. The 
second moment ^2 's perhaps a less familiar quantity to work with 
than the circular or line-of-sight velocities, but it has the distinct ad- 
vantage that it is not necessary to assume a particular 'splitting' be- 
tween ordered and random motions. We discuss its physical mean- 
ing in more detail in Section [5.5.3| 

The quantity b is often expressed as an anisotropy parameter 
I5~ = I — I /b, where j3; = corresponds to isotropy. To give an idea 
of the kinds of anisotropics observed in real galaxies, |Shapiro et al.| 
( |2003| l found flattened velocity ellipsoids with 0.35 < jSj < 0.75 in 
the discs of Sa and Sb galaxies. [Cappellari, et al.|(2007[ > and |Thomas| 
|et al.| ( [2009 ) found more circular ellipsoids for which fi^ < 0.4 in 
fast-rotator ellipticals and SOs. 

If the anisotropy j3, is free then the predicted kinematics for 
each galaxy are a function of three parameters: {M/L)x, Mdm and 
P~. We find, however, that for our galaxies, which are rotation dom- 
inated, the predicted kinematics along the slit are not very sensi- 
tive to the particular choice of j8;. This means that we are unable 
to place stringent constraints on anisotropy using these data, but it 
also means that we are not vulnerable to serious systematic errors 
due to our assumptions about anisotropy. 



2.4 Summary of methods and assumptions 

To sum up our methods section, we find the stellar component of 
each mass model by assuming axisymmetry and a constant stellar 
mass-to-light ratio. To each model we add a dark halo that follows 
a spherically symmetric NFW profile and assumes the correlation 
between halo concentration and halo mass defined by equation (jT2|. 
The total mass model is a function of two parameters, the stellar 
mass-to-light ratio {M/L)x and the dark halo mass within the virial 
radius Mdm- It is expressed as a sum of Gaussians. 

The mass model is used to calculate an estimate of the ob- 
served stellar kinematics. This is done by solving the Jeans equa- 
tions under the assumption of constant anisotropy in the meridional 
plane, yielding the second velocity moment, jJL2- This predicted 
quantity is then compared to the observed root mean square ve- 
locity Vnns = (v^ + CT^)'/^. The two parameters of the mass model 
{M/L)x and M^y^, and the velocity anisotropy j3j are then adjusted 
until the predicted stellar kinematics match the observations, plac- 
ing constraints on those parameters. 

Having determined the parameters for each galaxy, we also 
compute circular the velocities Vc of the stellar, dark and total mass 
distributions using the numerical techniques described in Cappel- 
|lari| ( |2002) . provides an intuitive measure of the mass enclosed 
as a function of radius, which is proportional to v^. The ratio of the 
squares of the stellar and dark circular velocities is therefore equal 
to the ratio of the stellar and dark mass enclosed. Vc will also be use- 
ful in a forthcoming paper on the TuUy-Fisher relation (Williams, 
Bureau & Cappellari, in preparation). 

All mass-to-light ratios in this work are given in solar units, 
i.e. the mass-to-light ratio at waveband X, (M/L)x, is in units of 
{M/L)x,Q. 



3 SAMPLE AND OBSERVATIONS 



We apply the methods described above to a sample of 28 edge- 
on disc galaxies selected by [Bureau & Freeman| ( [T999l l. 14 of the 
28 galaxies are classified as SOs and the remaining 14 are Sa~Sb 
galaxies (see Table[TJ. 

The galaxies in the sample were originally selected to investi- 
gate the nature of boxy and peanut-shaped (B/PS) bulges. Of the 28 
galaxies, 22 have a B/PS bulge and 6 form a control sample with 
spheroidal bulges. Galaxies with a B/PS bulges are edge-on disc 
systems with bulge isophotes above and below the centre of the 
disc that are either horizontal (boxy) or concave (peanut-shaped). 
B/PS bulges are thought to simply be bars viewed edge-on. Many 
A'-body simulations have shown that bars buckle and thicken soon 
after formation, resulting in an object that appears peanut-shaped 
when viewed side-on and boxy when viewed at intermediate an- 
gles ( e.g.|Combes & Sanders |198T||Combes et al.|1990||Raha et al.| 
[T99T]|Athanassoula & Misiri otis|2002^ . |Lutticke et al.| ( |2000| l found 
that 45 per cent of edge-on disc galaxies have a B/PS bulge, which 
is roughly consistent with the fraction of bars in face-on spirals. 
Moreover, several studies have identified the kinematic signatures 



of bars in systems with a B/PS bulge (e.g. [ Kuijken & Merrifield 
1995 ; Merrifield & Kuijken|1999l|Bureau & Freeman 1999 Chung 
|& Bureau 2004) and a recent study demonstrates the converse by 



observing kinematic evidence for a B/PS bulge in a face-on barred 
system jMendez-Abreu et al.|2008| l. The sample therefore contains 
galaxies whose central regions are probably non-axisymmetric. If 
there is a direct correspondence between B/PS bulges and bars, then 
the bar fraction in our sample (~ 75 per cent) is representative of 
that in all disc galaxies (~ 65 per cent, e.g. [Sheth et al.|2008| l. In 
fact, we will argue in Section|5]the these non-axisymmetries do not 
affect our results, despite our axisymmetric modelling techniques. 



To parametrize the projected light distribution we use the rel- 
atively deep, high-resolution A'n-band images presented by [Bureau] 
[eTar] i |2006) . In the course of this work we discovered that the 
photometric calibration of the images was incorrect (they were too 
faint by ~ 0.8 mag), so we used shallower 2MASS images to recal- 
ibrate them, transforming from Kn to A^^-band in the process. We 
describe this procedure in detail in Appendix[B] We also correct the 
2MASS images and therefore our models for the effects of fore- 
ground Galactic extinction using the dust maps of [Schlegel et al.| 
( |I998[ l. Because the apparent surface brightness distributions that 
we parametrize are at A^^-band, the mass-to-light ratios we measure 
are also at A^-band, and therefore denoted {M/L)k^. Throughout 
this work we adopt a value for the absolute magnitude of the Sun at 
A:5-band of Mks,q = 3.29 jBlanton & Roweis|2007^ . 



We compare the predicted kinematics of the mass models to 
major-axis long-slit stellar kinematics presented in [Chung & Bu-| 
|reau| ( |2004[ >. We use stellar kinematics rather than gas kinemat- 
ics because gas is not present and/or not extended in many cases. 
Where gas is present, it is strongly affected by non-circular motions 
and shocks in the irmer parts of many of the objects (see |Bureau &| 
I Athanassoulall 1 999\ |Athanassoula & Bureau||I999) . Line-of-sight 
velocity distributions were extracted using the Fourier Correlation 
Quotient algorithm ( |Bender|I990^ ; the v and a used are those of 
the best-fitting Gaussian. 
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Table 1. Galaxy sample 



Galaxy 


Type 


D 


^25 


Re 








B 




Me 






(Mpc) 


(arcsec) 


(arcsec) 






(mag) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 










B/PS bulges 












NGC 128 


SO pec 


59.2 


95.3 


18.1 


0.49 


2.58 


8.51 


12.46 


-25.35 


-21.40 


ESO 151-G004 


so 


106.9 


42.2 


14.1 


0.71 


2.13 


10.29 


14.73 


-24.85 


-20.41 


NGC 1381 


SAO 


16.8 


76.2 


20.1 


0.92 


3.47 


8.36 


12.35 


-22.76 


-18.77 


NGC 1596 


SAO 


15.3 


1 16.7 


21.7 


0.39 


2.11 


8.06 


11.94 


-22.86 


-18.98 


NGC 1886 


Sab 


24.5 


96.2 


31.1 


0.71 


2.20 


9.16 


12.22 


-22.79 


-19.73 


NGC 2310 


SO 


15.2 


135.6 


39.2 


0.64 


2.21 


8.43 


11.48 


-22.48 


-19.43 


ESO 311-G012 


so/a? 


14.5 


122.8 


29.7 


0.49 


2.04 


1.16 


10.89 


-23.05 


-19.92 


NGC 3203 


SA(r)0+ 


35.4 


96.6 


22.1 


0.64 


2.80 


8.86 


12.76 


-23.89 


-19.99 


NGC 3390 


Sb 


44.0 


93.8 


28.4 


0.92 


3.03 


8.39 


11.73 


-24.83 


-21.49 


NGC 4469 


SB(s)0/a? 


16.5 


143.3 


39.1 


0.60 


2.21 


8.09 


12.07 


-23.00 


-19.02 


NGC 4710 


SA(r)0+ 


16.5 


131.6 


44.0 


0.73 


2.19 


7.62 


11.69 


-23.47 


-19.40 


PGC 4493 1 


Sbc 


61.2 


91.4 


32.3 


0.66 


1.86 


9.60 


12.60 


-24.34 


-21.33 


ESO 443-G042 


Sb 


35.8 


84.6 


35.9 


0.86 


2.02 


9.37 


12.56 


-23.40 


-20.21 


NGC 5746 


SAB(rs)b? 


30.4 


217.3 


51.0 


0.67 


2.84 


6.86 


10.11 


-25.55 


-22.30 


IC 4767 


S pec 


53.3 


57.8 


21.4 


0.64 


1.72 


9.94 


13.91 


-23.69 


-19.72 


NGC 6722 


Sb 


70.9 


87.9 


22.4 


0.81 


3.19 


8.95 


12.28 


-25.30 


-21.97 


NGC 6771 


SA(r)0+ 


64.6 


98.7 


15.5 


0.52 


3.29 


8.97 


13.24 


-25.08 


-20.81 


ESO 185-G053 


SB pec 


66.9 


37.4 


12.3 


0.77 


2.34 


9.91 


14.04 


-24.22 


-20.09 


IC 4937 


Sb 


72.4 


51.5 


27.7 


0.87 


1.63 


9.70 


13.84 


-24.60 


-20.46 


ESO 597-G036 


SO" pec 


123.5 


43.3 


18.5 


0.87 


2.04 


10.10 


14.76 


-25.36 


-20.70 


IC 5096 


Sb 


47.0 


105.7 


27.6 


0.56 


2.14 


8.52 


12.15 


-24.84 


-21.21 


ESO 240-GOll 


Sb 


41.4 


167.2 


41.5 


0.73 


2.92 


8.43 


11.53 


-24.65 


-21.56 










Control sample 












NGC 1032 


so/a 


37.0 


106.0 


21.0 


0.58 


2.91 


8.40 


12.13 


-24.44 


-20.71 


NGC 3957 


SAO+ 


24.1 


100.7 


31.6 


0.65 


2.08 


8.62 


12.71 


-23.29 


-19.20 


NGC 4703 


Sb 


72.0 


50.9 


28.0 


1.33 


2.43 


8.92 


13.37 


-25.36 


-20.92 


NGC 5084 


SO 


24.2 


297.9 


27.0 


0.30 


3.26 


7.19 


10.07 


-24.73 


-21.85 


NGC 7123 


Sa 


55.3 


80.2 


17.1 


0.60 


2.83 


8.45 


12.90 


-25.26 


-20.81 


IC5176 


SAB(s)bc? 


24.8 


135.2 


29.6 


0.46 


2.10 


8.77 


12.02 


-23.20 


-19.95 



Notes. (1) Galaxy name. To ensure continuity with previous studies, the sample is subdivided into B/PS bulges and control galaxies and arranged in order of 



increasing right ascension. (2) Morphological type taken from|Jarvis| (T986) , |deSouza & Dos Aigos]{T987).|Shaw|fT987) and |Karachentsev et al.|(T993) . (3) 



Distance. The distance of NGC 1381 is taken from |Jensen et al. |^2003 and that of NGC 1596 from |Tonry et al.H200 lV NGC 4469 and NGC 4710 are members 
of the Virgo cluster so we adopt the cluster distance derived by Mei et al. (2007). For all other galaxies we use distances from the NASA/IPAC Extragalactic 
Database (NED) calculated assuming a WMAP5 cosmology I Komatsu et al. 2009 Hq = 70 km s~ * Mpc^' ) and a Virgocentic flow model I Mould et al. 2000 1. 
(4) Radius of the 25 mag arcsec"^ B-band isophote taken from HYPERLEDA ( Paturel et al. 2003 1. (5) Semi-major axis of the ellipse containing half the light 
at Ks-hanA, measured using the near-infrared photometry presented in |Bureau et al. (2006) and the method described in Section[4] (6) and (7) Radius of the last 
stellar kinematic data point used in this work, expressed as a fraction of R25 and R^. (8) Apparent magnitude at Kg-hanA derived from the MGE parametrization 
of the Kn-banA image jBureau et al. 12006} calibrated using the 2MASS Ks-hanA image jSkrutskie et al.|2006| connected for Galactic extinction. (9) Apparent 
magnitude at B-band corrected for internal and Galactic extinction taken from HyperLEDA. (10) and (1 1) Absolute magnitudes derived from the distance and 
apparent magnitudes in columns (3), (8) and (9). 



4 RESULTS 



4.1 Exploring parameter space 

Before presenting our best models, we give an overview of how the 
changes in {M/L)ks and Mqm affect the predicted second velocity 
moment jl2, and what kinds of constraints this allows us to place on 
these parameters. In general, increasing {M/L)k^ shifts the model 
/X2 to greater velocities at all radii. Increasing Mdm increases fl2 at 
large radii more than it does at small radii. Fig.[3]demonstrates this 
behaviour for NGC 1381. 

We define the best-fitting model for each galaxy to be the one 
with the combination of parameters (M/L)k^, Mdm and jSj that re- 
sults in a predicted ^2 that most closely matches the observed Vims- 
The figure of merit, is therefore defined as 



-Ai2[MDM,(M/L);f,,jS,]}2 



Av2 



(18) 



where Avrms is the error in Vims and the summation is over all kine- 
matic data points. 

We show contour plots of as a function of these pa- 
rameters for the complete sample in Fig. |4] For clarity, these 
plots are marginalized over the anisotropy P^. This means that the 
value of x^ used at each point in M-[)m~{M / L) Ks space is set to 
mm{x^[MoM,{M/L)Ks,l5z]) for < jS, < 0.75. The contours be- 
come horizontal for Mdm W^Mq, because very low mass haloes 
do not affect the dynamics of the galaxies. There are two classes of 
X^ plot: those where we can rule out = at at least the 3 (7 
confidence level (i.e. where the red contour in Fig. |4] fully encir- 
cles the minimum in x^, e.g. NGC 128 and NGC I38I) and those 
where we cannot (where the red contour forms an open horizontal 
'tongue', e.g. NGC 1886 and IC 4937). All galaxies therefore have 
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Figure 3. Predicted kinematics for several mass models of NGC 1381, demonstrating the effects of varying the model parameters Mum and {M/L)^^ Each plot 
shows the observed v,,„s (points with error bars) and predicted /i2 (solid red line). We also show the circular velocity of the luminous (dashed blue line) and 
dark (dotted blue line) components. The text at the bottom of each plot gives the parameters of the mass model in solar units at A'^-band. The best fitting mass 
model for this galaxy has {M/L)ks = 1.22, Mqm = lO'^'^' and ji, = 0.14. For this figure we assumed isotropy, i.e. /3j = 0.0. 



a best-fitting combination of mass model parameters and allow us 
to place upper and lower bounds on {M/L)ks and at least an upper 
limit on Mdm- Not all galaxies allow us to place a lower limit on 
Mdm- 



4.2 Best-fitting models 

The parameters of the best-fitting mass models are shown in Ta- 
ble |2] The MGE parametrizations of the luminous component of 
each mass model are shown in Table [aT| The complete results of 
the mass modelling and dynamical modelling procedures described 
in Section |2] are shown in Fig.[5](B/PS bulges) and Fig. |6] (control 
sample). In each panel: 

(i) The top plot shows the image and the luminous component of 
the mass model. The filled contours are the observed light distribu- 
tion while the solid lines are the isophotes of the best-fitting MGE 
model, projected and convolved with the PSF. Contours are sepa- 
rated by 0.5 mag arcsec^^. The lowest contour is 7 mag arcsec^^ 
below the highest. The horizontal axes of these images are regis- 
tered with those of the kinematic data below, which usually requires 



that the image be cropped. The full extent of the photometric data 
is presented in fig. 1 of |Bureau et al.| ( (2006| l. 

(ii) The bottom left plot shows the observed line-of-sight veloc- 
ity v (points with error bars) and the circular velocities of the to- 
tal mass model (solid blue line), the luminous component (dashed 
blue line) and the dark halo (dotted blue line) of the best-fitting 
mass model (see below). We remind the reader that the ratio of the 
squares of the dark and luminous circular velocities at a particular 
radius gives the approximate ratio of dark to luminous matter en- 
closed within a sphere of that radius (this statement is strictly true 
for spherical models only). 

(iii) The bottom right plot shows the observed root mean square 
velocity Vi-,^s (points with error bars), the second velocity moment 
jl2 of the best-fitting mass model with a dark halo (thick solid red 
line) and the second velocity moment of the best model without a 
dark halo (thin solid red line). 

We mark the radial axis of the bottom left plot with Ris /2 and 
^e, two distance scales used for spiral and elliptical galaxies, re- 
spectively. We place a mark at Ri^ /2 because i?25 is always beyond 
the extent of the stellar kinematic data. Showing both scales aids 
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Figure 4. contour plots for the complete sample showing fit quality as a function of dark halo mass and stellar mass-to-light ratio {M/L)k^ . The third 
parameter, the anisotropy has been marginalized out for this figure, as described in Section pTT] The red contours are the 3 a confidence levels. The inner 
black contours are at 1 fj and 2 fj. The cross in each plot shows the location of minimum x^- The upper horizontal axis shows the halo concentrations Cyji (Mqm) 
corresponding to Mdm. as defined by equation jl2^ . With the exception of ESO 240-GOll, we arrange galaxies in the same order in which we present them in 
Tablefl] The parameters of the best-fitting model for ESO 240-GOl 1 are outliers, so this galaxy is shown over a different range of {M/L)k^. 
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Figure 5. Complete results of the mass and dynamical modelling described in Section[2]for the sample of galaxies with a B/PS bulge. Briefly, the top plot of 
each panel shows the isophotes of the image (filled contours) and the MGE parametrization of this luminous component (solid black lines). The bottom left 
plot shows V, the observed line-of-sight velocity (points), v^, the total circular velocity (solid blue line), and the circular velocities of the luminous component 
(dashed blue line) and dark halo (dotted blue line). The thick notch on the S-axis shows R^, the X^-band effective radius. The thin notch shows S25/2, half 
the B-band optical radius. The bottom right plot shows Vi,„s, the observed root mean square line-of-sight velocity (points), 1I2, the best-fitting model second 
velocity moment (thick solid red line) and the best-fitting second moment with no dark halo (thin solid red line). See Section [4!2] for a more detailed explanation 
of these figures. 



comparison with previous studies of the influence of dark matter in 
both classes of galaxies. This is especiaOy interesting because our 
sample contains many SOs, which represent a transition between 
spirals and eUipticals. 

We note that the effective radius of very flat objects is partic- 
ularly sensitive to its definition. We define it here to be the semi- 
major axis of the ellipse containing half the integrated A'^-band 
light, which we determine directly from the photometry. We use 
the integrated light of the MGE model to determine the half-light 
radius rather than the more usual growth curve method. Using our 
MGE models, we tested that this definition yields a value similar 
to that which would obtained from an ideal observation with the 
galaxy oriented face-on. We found i?e(face-on) = 0.9i?e (edge-on) 
±0.1. For this definition of Rg and for our sample of edge-on spi- 
rals and SOs, R25 ~ 3.6i?e (see Fig.|7j. This empirical finding will 



be used in this paper, and is probably true to within approximately 
50 per cent for other edge-on objects, but it should obviously not be 
applied to less inclined spirals. 



An alternative classic definition of is to measure the radius 
of the circle containing half the light, independent of the apparent 
ellipticity l |Burstein et al.|1987| [Roman et al.|1991[|Jorgensen et ^ 
|1996^ . We measured this value for our edge-on galaxies and found 
that it was typically smaller than the semi-major axis value by a 
factor 2.0 ± 0.2. It is therefore rather poorly related to the 'true' 
face-on value and we prefer to integrate within ellipses. It is impor- 
tant, however, to keep our definition of Rf. in mind when comparing 
our results to other studies. 
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Figure 5. — continued 
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Figure 6. As Fig.|5]but for the control sample. 
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Table 2. Mass and dynamical modelling results 



Best-fitting models with halo Best-fitting models without halo 

Galaxy (M/L)jf; log(MDM/M-;,) ft ^red (M/L)j^^.nohalo ft.nohalo ;i;red.nohalo 



B/PS bulges 



NGC 128 


1 04+°°'' 

1.0t_Q Qg 


13 61+°"^ 


n 93+0.06 

U.ZJ q 


1.86 


1 9n+0.06 

^■^°-0.06 


17+0.08 

-0 09 


2.61 


ESO 151-G004 




12 S1+" " 




1.16 


1 59+0.11 

^••''-0.11 


(,21+0.15 
"■'^'-0.21 


1.38 


NGC 1381 


1 18+006 


13 21+"'3'5 
—0.40 




1.89 


1.48l!J-°l 


17+0.06 

-0.08 


2.85 


NGC 1596 


1 26+°-"' 

^^—0 14 


14 50+0.66 

—0 99 


14+"." 
—0.14 


1.17 


1.58l!J;|l 


12+"" 

—0 12 


1.80 


NGC 1886 


i.09;«i3 


12.53i';22 


4S+" 26 
o.tJ_o 43 


1.04 


1 C7+0.16 
^•-^'-0.15 


00+"'3'' 
"■""-0.00 


1.12 


NGC 2310 




13 26+" '* 

lO.zo 027 


o.ooliJig 


1.55 


1 49+0.07 
^•^'-0.06 


00+" "3 
"■""-0 00 


3.09 


ESO311-G012 


0.82;«'|« 


13.39;'j''^^ 


0.32l»|^ 


1.00 


99+0.07 

-0.07 


0.24lU-^2 


1.40 


NGC 3203 


0.77^«o^ 


13.17^!Ji« 




1.92 


1 17+0.03 

—0 03 


o.oolUSg 


3.94 


NGC 3390 




12.18+^:^5 


o.ooijj;;;^ 


1.33 


1.06lg'°l 


o^oolU;"^ 


1.63 


NGC 4469 


1 12+"" 

-0.11 


19 97+0.40 
-0.59 


05+"''' 
-0.05 


1.44 




00+" "' 
■^ -0.00 


1.81 


NGC 4710 


Q 7, +0.06 
-0.06 


,2 31+0.31 

-0.40 


00+" 
-0.00 


1.08 


95+".03 


00+" "2 

-0.00 


1.72 


PGC 44931 


1 04+°" 

^•^^-0.16 


19 04+0.36 


00+" ''' 
"•""-0.00 


1.75 


1 66+" '3 

-0.12 


00+" 
"■""-0.00 


2.79 


ESO 443-G042 


, 21+0.45 
^•^'-0.57 


1 1 gs+'.22 


47+".28 
'J-^'-0.47 


2.05 


^■°8_o.l8 


29+"'2'' 

"■"-0.29 


2.10 


NGC 5746 


, 9^+0.07 
^■^-'-0.05 


13 3S+".1« 


00+" "2 
"■""-0.00 


2.13 


1 64+""'' 
^■"^-0.03 


00+" °" 
"■""-0.00 


3.08 


IC 4767 


1 36+""* 
^-'"-0.16 


11 45+"'2 


24+"-"' 
"■"^^-0.23 


1.23 


1 51+"'3 
^■■^ -0.12 


l7+o.'7 

"■"-0.17 


1.39 


NGC 6722 


1 14+0.10 


1 2 79+".2'' 


n 14+0.09 


2.22 


1 47+0.06 
^•^'-0.06 


90+0.08 
"■'^"-0.11 


2.68 


NGC 6771 


1 19+0.12 
^•^^-0.12 


1 ^ on+O 68 


n 9Q+0 11 

"■^^-0.18 


2.05 


1 '11+0.10 
^•-'^-0.10 


n 19+0.17 
U.1Z_Q 12 


2.34 


ESO 185-G053 


1 n-j+oio 

^•"^-0.10 


1 on+0.81 
1^-^^-1.04 


0.34l»;?» 


1.73 


1 1 c+O.OS 

^•'-"-o.os 


"■zil!!:!^ 


1.93 


IC 4937 


1 07+"--^ 


J J 45+0.90 


0.18l!J;?^ 


1.35 




0.08lU;2« 


1.38 


ESO 597-G036 


1 01+0.41 


13.30+'J]^ 


17+0.26 
"■^'-0.37 


0.96 


1^87^oio 


00+0 40 
"■""-0.00 


1.44 


IC 5096 


1 n<;+0.10 


1 o ^T+O 50 

12.62^03, 


n rin+0.06 

"■OO-o.oo 


1.28 


1 on+0.05 

1 -2910.06 


n rin+0.05 

"■OO-o.oo 


1.70 


ESO240-G011 


2 41+0.14 
"'•^^-0.14 


10.561"'" 


00+" 
"■""-0.00 

Control sample 


1.76 


9 49+0.08 
^•^^-0.08 




1.76 


NGC 1032 


1 10+".18 

^•'"-0.12 


19 OS +0.63 
12.S5_05(, 


00+" -' 
"■""-0.00 


0.98 


1 34+0.16 
^•-'^-0.11 


0.08lg;i« 


1.42 


NGC 3957 


1 35+"''2 
^•-'-'-0.12 


■ 2 27+0.49 
'^•^'-0.63 


00+" -' 

"■""-0.00 


1.22 


1 65+""' 
^•"-"-0.06 


00+" ''' 
"■""-0.00 


1.63 


NGC 4703 


rv oo+O.lO 


1 2 40+0 36 
iz.tu 045 


09+" 23 
"■"^-0.09 


1.07 


1 06+" "** 


l2+"'2' 
u.iz u 12 


1.65 


NGC 5084 


0.89+^:^ 


13 03+" '''' 


00+" 
"■""-0.00 


1.18 


1 02+""*' 
^■"^-0.03 


00+" 
"■""-0.00 


1.68 


NGC 7123 


0.84^g:°? 


12.62+?;«i 


09+" '5 
"■"^-0.09 


1.18 


95+".05 

"■^-"-0.05 


06+"'"' 

"■""-0.06 


1.43 


IC 5176 




12.491'.;'' 


34+".23 
"■-'^-0.34 


1.38 


1 33+0.11 
'■-'•'-0.10 


06+" 2' 

"■"°-0.06 


1.44 



Notes. All values are computed assuming the distances in Table[T]and are presented in solar units at ^T^-band, (M/L)jfj EiTors are the formal fitting errors at 
the 3fT confidence level and neglect the uncertainties discussed in Section [53] If the lower error on halo mass Mdm is denoted with ellipsis, e.g. logMoM = 
12.531''"^, then for this galaxy there is no lower bound on Mdm because Mdm = cannot be ruled out at the 3 fJ level. In such cases the quoted value for Mdm 
is highly unconstrained up to some upper limit, and shoul d be used with caution. Xiei ™d ;i;ied, nohalo are the reduced x for the models with and without dark 
haloes, i.e. (;i;-/DOF) '/- where is defined in equation jisj and DOF is the number of degrees of freedom. 



5 DISCUSSION 

The form of the observed 

i'rm,s 't ftinction of radius is tistially 
a double-humped curve (with the inner hump reaching a smaller 
speed), but sometimes a monotonically rising curve reaching a 
plateau toward the edge of the disc (rather like a typical observed 
rotation curve). For example, the observed second moment rises 
monotonically in NGC 1886 and IC 5176, is almost constant in 
NGC 1032 and NGC 5084, falls before rising again in NGC 1381 
and NGC 2310, and rises then falls then rises again in NGC 3203 
and NGC 6771. The position, height, depth and number of bumps 
is different for each galaxy. The first thing to note, therefore, is that 
the models are able to reproduce this whole range of behaviours. 
This is demonstrated by Xied^ the square root of the reduced of 
the best-fitting models (see Table|2]l. The median x^^^ is 1.35 with 
an rms scatter of 0.40. 



These low figures of merit are neither trivial given the wide 
range of kinematics observed, nor necessarily expected given the 
simplicity of the models. Although the MGE technique allows for 
accurate parametrization of the A'^-band photometry, the depro- 
jected mass models are axisymmetric descriptions of objects that, 
in most cases, we have good reason to believe are in fact barred 
(see [Bureau & Freeman|p^99[ [Chung & Bur eau 2004). Our one- 
parameter description of the dark halo is model dependent, and our 
total model contains only three free parameters. Despite this, we 
are able to accurately reproduce the wide range of observed second 
velocity moment profiles as far as outermost data point, which is 
typically at .Rmax ~ 2-3 Rg or, equivalently, i?niax ~ 0.5-1 R25 (see 
Table [T}. Because of the lack of freedom that we have to vary the 
shape of the predicted kinematics, it seems that a great deal of in- 
formation about the kinematical structure of these early-type disc 
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Figure 7. The B-band optical radius R25 shown as a function of the Sj-band 
effective radius (as defined in Section [4!2) for our sample. The line is 
a linear fit to the data with the intersect constrained to be and the outlier 
NGC 5084 excluded. The best fit is R25 = 3.6Re. Symbols are as in Fig.[2] 
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galaxies is contained in their photometry alone. This result is con- 
sistent with that of Cappellari ( 2008 1: 

Given the extent to which our models reproduce the observed 
kinematics, we are justified to consider that they correspond in some 
way to the intrinsic properties of the galaxies. We defer our discus- 
sion of the possible systematics introduced by our assumptions to 
Section l533] 

The following discussion of the model parameters is separated 
into two parts. In the first, we discuss the properties of the best- 
fitting mass models that include a dark halo. In the second, we 
discuss those without. We present the three parameters, {M/L)xs, 
Mdm and (67, of the models with dark haloes and the two parame- 
ters, {M/L)iCs and jS^, of the models without in Table|2]and Fig.js] 
We precede this discussion by noting than none of them appears to 
correlate with the absolute magnitude Mk^, the stellar mass or the 
Hubble type. 



5.1 Models with dark haloes: mass-to-light ratios 

We find a median A'j-band mass-to-light ratio {M/L)ks = 1 -09 with 
an rms scatter of 0.31. As we noted in the introduction, the mass- 
to-light ratios are themselves interesting because there are relatively 
few dynamical measurements of stellar K-hand mass-to-light ratios 
and they are important for normalizing and constraining stellar pop- 
ulation models in the near-infrared. We therefore compare our val- 
ues to other available dynamical and stellar population-based mea- 
sures. Of course there have been many measurements of (M/L) ra- 
tio at wavebands other than Kg that can be transformed to using 
a global colour. To minimize the uncertainties associated with such 
transformations, we restrict our discussion to direct determinations 
at H, K or, almost equivalently, A'^. K — Kg < 0.03 mag (the exact 
value depends on which A'-band filter is used). For a given stel- 
lar mass, this is equivalent to a difference between {M/L)ks and 
{M/L)k of less than 3 per cent. We neglect this difference in the 
discussion that follows. 




Figure 8. Distributions of the parameters of the best-fitting models with 
a dark halo (solid black) and without a dark halo (hatched white) in Ta- 
ble |2] We do not show halo masses where there is no lower mass limit. 
Corresponding values of concentration c^\j[M\}y{] are shown on the upper 
horizontal axis of the halo mass histogram. 



5.1.1 Comparison to previous dynamical mea.sures of stellar 
{M/L)k 



[Devereux et al.| jl987[ l made among the earliest determinations 
of near-infrared mass-to-light ratios using nuclear dispersions and 
1.65 lim (i.e. //-band) photometry. For 72 bright spiral, SO and el- 
liptical galaxies, they found a mean {M/L)h = 0.94 (where we have 
updated the value they adopted for Mq^h with a more recent deter- 
mination). For a constant characteristic H — Kg color of 0.28 mag 
I Jarrett et al.|2003| , this implies {M/L)ks ~ 0.7. |Kranz et al.||2003 1 
determined maximal disc {M/Vjg ratios by comparing the results 
of hydrodynamical simulations to the observed spiral patterns in 
spiral galaxies. For a sample of five high surface brightness late- 
type spiral galaxies, they find a mean (M/L)k ~ 0.6. In the case 
of the small elliptical galaxy Cen A (C7e ~ 140 km s^'), there are 



two independent dynamical measurements of {M/L)k |Silge et al 
2005 . Cappellari et al. 2009 ), both of which are ^ 0.7. We therefore 
find that our direct dynamical measures of {M/L)x^, which we be- 
lieve accurately account for dark matter, are somewhat larger than 
most previous dynamical estimates at around the 1-2 CT significance 
level. Our method depends on an assumed single-parameter NFW 
dark halo, rather than the two-parameter isothermal halo used by 
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Figure 9. Stellar (M/VjK^ against B — Ks colour. The data points use the 
A's-band mass-to-light ratios of the stellar components of the best-fitting 
mass-models including a dark halo, i.e. column (2) of Table |2] and B-band 
magnitudes, coiTected for Galactic and internal extinction, are taken from 
HyperLEDA. Kg-hand magnitudes are derived from our MGE models which 
incorporate a correction for Galactic extinction. Internal extinction at K- 
band is neglected. A median error bar that does not include the uncertainties 
introduced by the extinction corrections is shown in the upper left corner 
(see text). Symbols are as in Fig. [2] The thin straight line is the prediction 
of the spectrophotometric galaxy evolution models of Bell et al. 1 2003i. The 
thick straight line is a fit to our data constrained to the slope of the dashed 
line. The evolutionary tracks for a range of metallicities are from the models 
of |Maraston|p005) and run from 3 Gyr (low mass-to-light ratio) to 15 Gyr 
(high mass-to-light ratio). Metallicity increases from left to right for these 
tracks. 



|Kranz et al.H2003] l. It is however somewhat more direct and avoids 
the maximal disc assumption. 



5.1.2 Comparison to previous stellar population measures of 
{M/L)k 

If our modelling assumptions are justified, then because mass is 
correctly shared between the luminous and dark components, the 
{M/L)ks ratios that we measure should agree with those predicted 
by well-calibrated stellar population models. A forthcoming paper 
will present an extensive stellar population analysis of the galaxies 
in our sample using absorption line strengths and a more complete 
comparison to the predictions of models. For now though, we per- 
form a preliminary check by comparing our results to the predic- 
tions of two models. In Fig.|9]we show the dynamically determined 
stellar mass-to-light ratios of our sample galaxies as a function their 
B — Ks colour. We also show the predictions of two independent 
stellar population models. 

We use the S-band magnitudes given in HYPERLEDA, which 
are corrected for both Galactic and internal extinction. We adopt 
the uncertainties in the uncorrected 6-band magnitude, which ne- 
glect random and systematic errors in the extinction corrections 
and are therefore lower limits. The Galactic correction is derived 
from the dust maps of |Schlegel eTaLj l fTggS} . The median Galactic 
correction at B is 0.2 mag and the uncertainties are probably rather 
smaller than this. ESO 311-G012 is close to the equatorial plane of 
the Galaxy (Galactic latitude b = —8.0°) so the Galactic correction 
for that galaxy is very large (1.7 mag) and is derived from a re- 
gion of the |Schlegel et al.| ( |1998[ l maps where extinction corrections 
are unreliable (b < 5°). The B-band internal extinction coiTection 



is based on the statistical method of |Bottinelli et al.|jl995| l, which 
parametrizes the correction as a function of axial ratio (and there- 
fore inclination) and bulge-to-disc ratio (and therefore morpholog- 
ical type). In exactly edge-on spiral galaxies, the coiTections used 
are both large (0.8-1.5 mag) and somewhat dubious. The A'x-band 
magnitudes include a small correction for Galactic extinction (the 
median correction is 0.02 mag) and neglect internal extinction. The 
effect on the parameters of the model of neglecting internal extinc- 
tion in the A'^-band is discussed in Section [5. 5. 2[ but as far as the 
B — Ks colour is concerned this should introduce a systematic but 
small error (perhaps 0. 1 mag in the spirals and less in the SOs). We 
therefore conclude that the uncertainties in the B — K colours are 
dominated by large uncertainties in the corrections to the i?-band 
magnitude, especially the internal correction. 

The first prediction we compare our results to is the colour- 
(M/L) relation of Bell et al. ( 2003) . This relation is derived from 
spectrophotometric galaxy evolution models that use the PEGASE 
stellar population model (see |Fioc & Rocca-Volmerange|I99'7] and 
[Bell et al.|2003^ with a Salpeter initial mass function (IMF) mod- 
ified by globally scaling down the stellar mass by a factor of 0.7. 
|Bell et aL] l |2003^ do not present relations involving near-infrared 
colours, so we adopt a characteristic R — Ks colour of 2.2 mag to 
derive the line shown in Fig. |9] from their {B — R)-{M /L)ks 



tion, y ielding log(M/L)j^, = 0.138 -0.568(g-.g5). |Bell & de Jong 
( |200I| l do use near-infrared colours, but their relations are based 
on a small selection of ages and metallicities rather than the distri- 
bution observed in the local universe. This leads them to derive a 
rather flatter relation than the ones shown in |Bell et al.| ( |2003y 

We show in Fig.|9]a further comparison to an independent set 
of stellar population models from [Maraston| (IQQS) that predict A"- 
band mass-to-light ratios for stellar populations with ages less than 
15 Gyr and metallicities — 1 .33 < [Z/H] < 0.67. The models include 
a semi-empirical treatment of the thermally pulsing asymptotic gi- 
ant branch, to which the near-infrared luminosity is particularly sen- 
sitive. We show models assuming a Kroupa IMF. These tracks can 
be transformed to predictions for the scaled down Salpeter IMF of 
[Bell et aL| ( [2003[ > and the normal Salpeter by adding 0.09 dex and 
0.24 dex respectively to the logarithm of the (M/L) ratio. In this 
case we show evolutionary tracks for a range of metallicities rather 
than a statistical characterisations of local galaxies. 

As the figure demonstrates, the predictions of the [Bell et aL| 
( |2003| l galaxy evolution models are below our dynamical determi- 
nations in all but four cases. A power law fit to our models, con- 
strained to the same slope as the |Bell et aL| |2003| line and from 
which the outlier ESO 240-GIl is excluded, implies a systematic 
offset of Alog(M/L);fj = 0.06 dex. At least part of this small offset 
could be due to the introduction of a systematic error in the approx- 
imate transformation of the [Bell et al.| p003') io B — Ks. Our mea- 
surements are systematically toward the upper end, if not outside, 
the range of values predicted by the |Maraston| ( |2005) models. If the 
stellar population models are correct, this implies systematically old 
stellar populations, > 10 Gyr. These differences are intriguing and 
could suggest problems in the stellar evolution models or the IMFs 
adopted, but we note four reasons which lead us to argue that our 
findings are broadly consistent with their predictions. Firstly, there 
is a random scatter of O.I- 0.2 dex in the |Bell et al.H20()3) relation 
(larger at the blue end). Secondly, there may be systematic errors 
in the (B — Ks) colours adopted for our sample of edge-on galaxies 
due to internal extinction. These are difficult to quantify. Thirdly, 
in the case of the |Marastori| p005[ > tracks, the mass-to-light ratios 
are for single stellar population models, while many galaxies, es- 
pecially spirals, are likely to have composite stellar populations. 
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Finally, our sample of 28 bright early-type disc galaxies is not nec- 
essarily representative of local galaxies. We will perform a more 
detailed comparison on the models and data in a future work. 

We end this discussion of the (M/L)ks ratios measured by not- 
ing the existence of a significant outlier galaxy, ESO 240-GOll, 
which has {M/L)ks = 2.41 ±0.14. This is more than four standard 
deviations above the median (M/L)^^ of the sample. ESO 240G11 
has the latest Hubble type in the sample, and later galaxies are gen- 
erally thought to be more dark matter-dominated. However, our 
modelling method aims to correctly account for dark matter so 
{M/L)ks should be the mass-to-light ratio of the stellar component, 
not the total matter distribution. This particular galaxy therefore re- 
mains puzzling. 



5.2 Models with dark haloes: dark halo masses 

For those galaxies with constrained halo masses, we find a me- 
dian = lO'^'^^Mpj with an rms scatter of 0.7 dex. Using the 
concentration-halo mass correlation of |Macci6 et al.|p008| >, which 
our haloes are constrained to lie on, this corresponds to Cvi,- = 7.9 
with an rms scatter of 1.2. These results should be treated with cau- 
tion, however, because most of the dark mass lies beyond radii at 
which our models are constrained by either photometry or kine- 
matics. Mdm is therefore strongly dependent on the assumptions 
described in Section lT2l 

With this strong caveat in mind, we note, however, that the 
virial masses Myjj = \.2My)m of the models are at least consistent 
with the predictions of semi-analytic and halo occupation distribu- 
tion models with respect to their stellar mass. We present a compar- 
ison in Fig. [To] Croton et al.|p006) and |Somerville et al.| ( |2g08) 
use semi-analytic models of galaxy formation that, among other 
things, attempt to predict the relationship between the stellar mass 
and dark halo mass of a galaxy by modelling relevant physical pro- 
cesses such as the growth of structure, cooling, star formation and 
feedback. [Wang et al.| j2006[ > and |Moster et al.|p009) use a halo 
occupation distribution approach to populate a halo catalogue with 
stellar mass. The halo catalogue is taken from |Springel et af]j2005^ . 
It is populated to reproduce the local observed stellar mass function 
and clustering function. 

Our models do not agree with the predictions in a few cases, 
but there is no evidence of any systematic offset. The stellar and 
virial masses of our mass models do not seem to be correlated. This 
is not true of the predictions, especially below = 1O'^M0, 
which we unfortunately do not probe. The comparison is there- 
fore limited by both a lack of mass range in our sample, and the 
significant uncertainties in our halo mass determinations. We can 
only conclude that our galaxies have stellar mass to halo mass ra- 
tios consistent with the range predicted by the theoretical models 
for these halo masses. However, the assumptions and observational 
constraints used to make their predictions are quite independent 
from ours, so even this relatively weak statement is not trivial. Ig- 
noring the predictions, we note that the spiral galaxies in our sample 
all have large M*/Mvii- ratios relative to the average SO. 

The dark matter enclosed within a sphere of radius r is given 

by 



12 F 



AfCvir) ■ 



(19) 



From this and the ratios of the circular velocities of the dark and 
stellar components, we can calculate the dark-to-total mass fraction 
as a function of radius. We find a mean X^^lRg) of 15 per 
cent (10 per cent rms scatter) and a mean Xdm(^25) of 49 per cent 
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Figure 10. Comparison of the stellar and halo masses of our mass models 
(data points) to the independent predicti ons of [Wang et al.H2006) , [Croton| 
|etarH2006t , |Somerville etaflpOSg) and |Moster et al.H2009{ . Median error 
bars are shown in black in the upper left comer. Symbols are as in Fig. [2] 



(15 per cent scatter). As we mentioned in the introduction, a max- 
imal disc is strictly defined as one in which the stellar component 
contributes 75-95 per cent of the rotational velocity at 2.2R^[sc^ i-S- 
XDM(2.2.Rdisc) < 45 per cent jSackett|1997) (LlR^i^c is the radius 
at which the rotation curve of a thin disc reaches its maximum). The 
disc scale length of edge-on, barred systems are uncertain ( |Bureau| 
[et al.|2006| >, but our estimates of -Kdisc suggest that X£)y[{2.2R^j^^) 
is 27 per cent (11 per cent scatter). The best fitting mass model is 
sub-maximal in only two cases [NGC 2310, XDM(2.2i?(jisc) = 46 
per cent and PGC 44931, XDM(2.2.Rdisc) = 50 per cent]. 

We can also determine the radius at which dark and stel- 
lar matter contribute equally to the circular velocity and mass en- 
closed within a sphere. We do this in units of Rg and i?25 and find 
^(^DM = 0.5)/^e = 4.1 (with rms scatter of 1.3) and R{Xum = 
0.5)/iJ25 = 1.0 (with rms scatter of 0.5). These results should 
be treated with a little caution, however, because in several cases 
^(^DM = 0.5) > Rm-dx, the last kinematic data point. Beyond the 
kinematics, the models are obviously particularly model dependent. 

As discussed in Section 5.5.1 Xqm and R{X£)m = 0.5) are al- 



most unaffected by distance or photometric calibration errors. Fi- 
nally, we note that there is no evidence of any significant difference 
between the dark matter content of the spirals and SOs in our sam- 
ple, at least at the radii within which the models are constrained by 
the kinematic data. 



5.3 Models with dark haloes: anisotropies 

The best-fitting anisotropies for our inodels are in the range 0.0 < 
P- < 0.5. As we discussed in Section [23] however, the constraints 
on these measurements are weak because the galaxies in our sam- 
ple are rotation-dominated and we do not have integral field data. 
Moreover, we assume a constant anisotropy throughout the galaxy, 
which is of course a significant simplification, since many of our 
galaxies have a bulge, bar and disc. With these caveats in mind, the 
typical anisotropy we find is broadly consistent with the average of 
previous observations of ellipticals, that are expected to have similar 
dynamics to the bulges (j8; < 0.4, Cappellari, et a l.|2007[|Thomas| 
[et al.|2009| l and discs (0.35 < jS^ < 0.75, , Shapiro et al.|2003[ >. 
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5.4 Neglecting dark matter 

As discussed in tiie introduction, previous studies have shown that it 
is usually possible to construct mass models of bright spiral galaxies 
that reproduce the general form of their rotation curves within i?25 
with no or only a sub-dominant contribution from dark matter. It is 
therefore not surprising that, in some cases, the improvement in fit 
quality when a dark halo added is not overwhelming. 

How many galaxies can we fit without a dark halo? Many of 
the second velocity moments of the best mass models without a dark 
halo are, by eye at least, satisfactory fits to the observed kinematics 
(compare the thin red lines to the data points in the lower right plot 
of each panel in Figures|5]and|6]l. In the formal sense, in four cases 
(NGC 1886, IC 4937, ESQ 240-GOll, and IC 5176) the removal of 
the dark halo does not improve the fit at the 3 CT level. 

We therefore now discuss an Occam's Razor argument which 
states that, since we can fit some of these galaxies well without dark 
matter, dark matter is not a significant component by mass within 
the radii probed by our kinematic data (0.5-1 i?25 or 2-3 i?e) in at 
least those cases. We reject this argument because the parameters 
of our best fitting mass models with dark matter are in line with 
independent expectations, whereas those for the haloless models 
are not. For example, the mass-to-light ratios we measure for mass 
models without haloes are inconsistent with stellar population re- 
sults. We find the median (M / L)i^^ j^q-^^\q = 1. 25 (M/L)jf J, where 
(M/L)/fj nohalo is the stellar mass to light ratio measured using a 
mass model without a dark halo. The (M/L)/^^ nohalo measure 
is therefore significantly higher than the predictions the stellar pop- 
ulation models of |Bell et al.| j2003| and Maraston| j2005| . Either 
the stellar population models and their assumptions are flawed, or 
the stellar mass models without dark matter have too much mass. 
The halo sizes (and concentrations) of the models with dark matter 
also match the expectations of galaxy formation. Of course, adding 
another free parameter to a model is always going to improve the 
quality of the fit, but if the true halo of a galaxy is not significant 
(or not present), the agreement between the halo masses that most 
improves the fits and the predictions and the expectations of ACDM 
would be a striking coincidence. 

Indeed, the fact that the (M/L) ratio of a mass model includ- 
ing only stars systematically exceeds that of a mass model account- 
ing for dark matter or a stellar population estimate can be used to 
estimate how ' wrong' the former is. The median (^/^)^^,nohalo 
must be decreased by 20 percent [A log (M/L) ~ —0.10] to match 
(M/L)/j-j, the dynamically determined stellar mass-to-light ratio 
with dark matter. This is similar to the finding of [Cappellari et aL] 
1 2006[ l, that the median /-band dynamical (M/L) ratio neglect- 
ing dark matter must be decreased by 30 per cent [A log (M/L)/ ~ 
—0.15] to match the stellar population (M/L) ratio. 

5.5 Uncertainties and assumptions 

5.5.1 Random errors 

The errors presented in Table|2]are 3 CT formal fitting errors. These 
neglect random errors due to distance uncertainties and photomet- 
ric calibration. As shown in Table [T] we adopted surface bright- 
ness fluctuation distance estimates in two cases (NGC 1381 and 
NGC 1596) and the Virgo cluster distance in two others (NGC 4469 
and NGC 4710). The uncertainties on these estimates is likely to 
be ~ 1 Mpc, limited either by the surface brightness fluctuation 
method or the physical size of the Virgo cluster. All other distance 
estimates are redshift-based estimates from NED assuming a Virgo- 
centric flow model ( Mould et al..2000) . The uncertainties on these 



distances are much larger, probably ~ 20 per cent. Because of the 
calibration process described in Appendix |B] the error on the cali- 
bration of the A^^-band images we use is of order that of the 2MASS 
survey (Skr utskie et al.|2006j l, which is negligible compared to the 
distance uncertainties. 

Random photometric and distance errors propagate linearly 
into the uncertainty in total dynamical mass of the model. In the 
presence of a dark halo these errors propagate into the stellar and 
dark parameters of the mass model in a complex and non-linear 
way. However, we find empirically that a given fractional error in 
photometric calibration or distance propagates into the same ran- 
dom error in (M/L)^^ and Mdm to within a factor of 2 or so. In 
fact, the formal fitting errors overwhelm these observational uncer- 
tainties. This can be seen from Fig. [4] where the 3 (7 confidence 
region allows us to constrain Mqm within ~ 1-2 dex only, a much 
bigger uncertainty than the distance and photometry errors. 

Despite the non-linear way in which photometric distance er- 
rors propagate into the stellar and dark masses of the mass models, 
we note that these errors almost cancel out in XDf^{R), the dark- 
to-total mass fraction at radius R. For example, plausible distance 
errors of 20 per cent result in changes in X£,MiRe), -?''DM(2.2/?dise) 
and Xqm {R25 ) of less than 2 per cent. 



5.5.2 Systematic errors due to dust 

We now discuss two sources of systematic error, both of which are 
due to dust. Firstly, if the internal extinction at A^^-band in these 
galaxies is significant, then the current modelled mass-to-light are 
overestimates. If that absorption varies significantly across the pro- 
jected galaxies, then this would also affect the halo masses deter- 
mined. We believe that this effect leads to < 10 per cent reduction 
in the total light detected at A^^-band, and probably much less. Ab- 
sorption at A'5-band it typically smaller than at optical wavelengths 
by an order of magnitude. This is reflected in detailed examination 
of the R-hwA ( [Bureau & Freeman] 1 999[ > and A'^-band (Burea u et al.| 
|2006^ images of the sample. Dust lanes that are prominent in the op- 
tical cannot usually be detected at A^j-band and half of the sample 
are SOs, which do not have optical dust lanes at all, so are likely to 
be totally unaffected by dust at near-infrared wavelengths. 

The second potential systematic effect is due to the position- 
ing of the slit used to measure the stellar kinematics. As described 
in [Bureau & FreemanH1999^ and [Chung & Bureau[ ( f2004[ l, in the 
dustiest, exactly edge-on systems the long-slit had to be moved a 
little away from the major axis. If the kinematics at this distance 
from the plane are significantly different from that along the ma- 
jor axis, then it would likely affect both the (M/L) ratio and halo 
mass determinations. The most likely result would be a systematic 
reduction of both parameters, since both rotational velocity and dis- 
persion typically fall with height above the disc. In actual fact, how- 
ever, the effect may be smaller than this in the B/PS sample. B/PS 
bulges are believed to rotate cylindrically, i.e. their kinematics do 



not vary with height ([Jarvis[[I987 


Shaw et al.[[1993|[Shaw[[1993| 


[Fisher et al.[I994|[D'Onofrio et al. 


199911. If this is true then this ef- 



feet should be small within the B/PS bulges. More importantly, the 
sample was constructed to avoid placing the slit above the plane. 
The dustiest galaxies in the sample (e.g. NGC 5746) are in fact in- 
clined at slightly less than 90°. It was therefore almost always pos- 
sible for the centre of the galaxy to lie within the slit, without dust 
affecting observations much. 
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5.5.3 Systematic errors due to model assumptions 

Finally we discuss the systematic errors introduced by our mod- 
elling technique, which assumes that the stellar mass-to-light ratio 
and anisotropy are constant across each galaxy, that the halo is of 
the form described in Section [2^ and, most importantly that the 
galaxy is axisymmetric. 

In the case of both {M/L)ks and jS^, the best-fitting param- 
eters may be thought of as global averages for each galaxy. We 
are, however, justified in assuming a constant value in both cases. 
The stellar population mass-to-light ratio at A'^-band varies rela- 
tively little between galaxies ^Bell et aLl|2003[ l, and rather less 
within them. As discussed in Section |2.3| and Section [53| our results 
are largely insensitive to the anisotropy chosen because our galax- 
ies are rotation-dominated. Choosing a more sophisticated, varying 
anisotropy would not change our conclusions. 

It is one of the goals of this study to offer constraints on halo 
masses and dark matter fractions by assuming a particular, theoret- 
ically motivated, single-parameter halo model. We are effectively 
asking how large haloes are, not what shape they have. If real haloes 
do not follow the NFW profile then our estimates of its virial mass 
are unlikely to be correct. We are encouraged, however, by the good 
correspondence between the stellar masses of the mass models in- 
cluding a halo and the predictions of stellar population models. This 
is strong circumstantial evidence that our method correctly appor- 
tions mass to the stellar and dark components, despite (or perhaps 
because of) the strong assumptions we make about the halo. 

Conclusions drawn from these axisymmetric models may be 
flawed if the assumption of axisymmetry is significantly violated. In 
line with local disc galaxy demographics, a majority of the galaxies 
in our sample are thought to host bars (see Section [3](. The possi- 
ble consequences of the application of axisymmetric modelling to 
the present sample (or indeed any representative local sample) are 
therefore a concern. Indeed, there are hints of our modelling break- 
ing down, presumably due to a bar, in the inner ~ 10 arcsec of some 
galaxies, where the model second velocity moment overpredicts 
the observations (e.g. NGC 1381, ESO 240-GOll and NGC 3957). 
None of the models underpredicts the data in this region. This can 
easily be explained by the sample selection. The B/PS bulges dom- 
inating our sample are thought to be bars viewed side-on. As such, 
the orbits supporting the bars are elongated perpendicular to the 
line-of-sight, resulting in observed velocities smaller than the cir- 
cular velocities of the equivalent axisymmetric mass distribution. 
End-on bars, which lead to the opposite effect, are systematically 
excluded from our sample as they appear round. This slight but 
systematic failure of our models in the central regions is therefore 
consistent with previous studies of this sample ( Bureau & Freeman] 



T999'. 'Chung & Bureau"2004'; 'Bureau et al.'2006') and others (e.g. 
Kuijken & Merrifield 1995 , Merrifield & Kuijken 1999) supporting 



a close relationship between B/PS bulges and bars. 

This, however, does not mean that our axisymmetric mass 
models are not appropriate. Firstly, although a bar distorts the sur- 
face brightness of a galaxy, the potential and the shapes of the orbits 
themselves are significantly less distorted (see, e.g., figure 2 (a) of 
[Bureau & Athanassoula|1999[ which demonstrates for a particular 
mass model that the radial peak in the axial ratio of the orbits sup- 
porting the bar is quite narrow). It is these orbital shapes that the ob- 
served stellar kinematics depend upon. This probably explains why 
the region in which the over-prediction occurs is much smaller than 
the B/PS bulge. The effect of the bar, at least as far as our method 
is concerned, is restricted to only the inner few arcseconds. 

Secondly, there is no systematic difference between the fit 



quality of haloless models in the inner regions of the galaxies with 
B/PS bulges and the control sample of six galaxies with spheroidal 
bulges. Unfortunately, we are limited by the rather size of the con- 
trol sample. Demonstrating a statistically significant correlation be- 
tween /l2 being too large in the central few arcseconds and more 
reliable bar diagnostics (e.g. [Bureau & Freemml|I999[ [Chung &| 
[Bureau|2004[ l is therefore unlikely. 

Nevertheless, it seems that the pressing question is not so much 
why the stellar kinematics of so many galaxies deviate slightly but 
systematically in the central regions, but rather why Jeans mod- 
elling of axisymmetric mass distributions reproduces the kinemat- 
ics of barred galaxies so well. The high quality of the kinematic fits 
perhaps suggests that the jl2 profiles primarily trace the enclosed 
masses as a function of galactocentric radius and are fairly (but not 
entirely) insensitive to the details of the dynamics. Because of this, 
axisymmetric mass models that reproduce the 'bulge-plateau-disc' 
surface brightness profile often observed in this sample (see Fig.[T] 
and jBureau et al.|2l}06^ can also reproduce the second velocity mo- 
ment of intrinsically non-axisymmetric galaxies with the same ra- 
dial profile. Irrespective of the importance of pressure support in 
our sample galaxies, /X2 therefore appears to be analogous to the 
circular velocity in purely rotationally-supported systems. 

We further note that there is no physically motivated evolu- 
tionary scenario which would systematically lead to axisymmet- 
ric galaxies with such a radial mass distribution, but such surface 
brightness plateaus develop naturally in barred disc galaxies due to 
the angular momentum and mass exchanges mediated by the bar 
(see, e.g., [Bureau & Athanassoula|2005| l. In effect, observations of 
jl2 circularize the intrinsic mass distribution so jlj is simply too 
coarse a measurement to properly constrain the internal dynam- 
ics of galaxies. To do so requires knowledge of the full shape of 
the line-of-sight velocity distribution. This approach was used by 
[Chung & Bureau[p004) for this sample. 

In summary, we argue that the overall overwhelming quality 
of the fits for galaxies with and without B/PS bulges provides con- 
fidence in the reliability of the derived masses and mass-to-light 
ratios, despite the likely non-axisymmetry of the inner regions of 
many of the galaxies. Our results, however, are not inconsistent with 
observations and simulations that demonstrate that these galaxies 
are barred. The application of the JAM technique to model galaxies 
of known non-axisymmetric morphologies is of course the best way 
to confirm definitively the validity of this argument. 



6 CONCLUSIONS AND FUTURE WORK 

We presented mass models for a large sample of spiral and SO galax- 
ies. These models allowed us to constrain the stellar and dark mat- 
ter content of the sample galaxies. For each galaxy, the stellar mass 
distribution was derived from near-infrared photometry under the 
assumptions of axisymmetry and a constant stellar mass-to-light 
ratio. We added an NFW dark halo and assumed a correlation be- 
tween concentration and virial mass. We solved the Jeans equations 
for the corresponding potential under the assumption of constant 
anisotropy in the meridional plane. By comparing the predicted sec- 
ond velocity moment to observed long-slit stellar kinematics, we 
determined the best-fitting parameters of the mass models. In some 
galaxies the observed second velocity moment rises monotonically, 
in others it plateaus at small radii, and in others it falls significantly 
before rising again. Despite this wide range of observed behaviours, 
our simple models, with only three free parameters (stellar mass-to- 
light ratio, dark halo mass and anisotropy), are able to reproduce the 
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observed kinematics very well. The observed kinematics typically 
extend to 2-3 or, equivalently, 0.5-1 /?25, and 1 < Xred ^ 2 for 
all galaxies. 

For our sample of 14 spirals and 14 SOs, we find a median 
{M/L)ks of 1.09 {M/L)ks,q with rms scatter 0.31. Our values are 
roughly consistent with the small number of previous independent 
determinations using different dynamical methods. The Ks-h&nA 
mass-to-light ratios that we measure are unique, however, because 
of the size of the sample and the way we have attempted to cor- 
rectly apportion mass between the stellar and dark components of 
the galaxies, without resorting to either a maximal disc assumption 
or results from stellar population models. Because they do not de- 
pend on stellar population models, they can be used to attempt to 
constrain the normalization and IMF of such models in the near- 
infrared. 

We also performed preliminary comparisons of our dynamical 
(M/L) ratios to the predictions of two stellar population models: the 
color-(M/L) relations of|Bell et aL](|2003^, which use the PEGASE 



stellar population models, and evolutionary tracks in color-(M/L) 
space for a range of metallicities from |Maraston| pOOS V The 'Bell' 
|et al.| l|2003 ) prediction is offset from our models by a small but sys- 
tematic amount and some of our galaxies are just outside the range 
expected by the |Maraston| ( |2003) l models. These differences could 
be due to systematic errors introduced in our comparison, but may 
also hint at problems with the stellar population models or the IMFs 
they assume. In a future work we will extend this comparison by 
determining absorption line strength indices for the present sample, 
allowing us to directly compare dynamical and stellar population 
estimates of {M/L) for individual galaxies. 

Once accurately known, ATj-band mass-to-light ratios are par- 
ticularly useful for constraining the stellar mass budget of the uni- 
verse. Firstly, (or K) is the waveband at which the effects of 
dust on observations of light from stars is minimized. Shorter wave- 
lengths are subject to absorption and longer wavelengths have con- 
tributions from hot dust seen in emission. Secondly, A^-band is dom- 
inated by light from the sub-solar mass main sequence stars that 
dominate the total mass budget (and trace the smooth potential) of 
galaxies. It therefore comes as no surprise that, as found dynam- 
ically in this work and using stellar population models (e.g. |Bell| 
et al.||2003| , {M/L)k^ is a relatively constant quantity in the lo- 



cal universe compared to the {M/L) ratio at 6-band. We both find 
that {M/L)ks varies by < 0.3 dex (a factor of 2) across our sam- 
ples. This is much less than the variations observed at shorter wave- 
lengths. The A'j-band luminosity function can therefore be used as 
a reliable proxy (subject to the small variation in (M/L)^^) for the 
stellar mass density function (e.g. [Bell et al.|2003| [Amouts et al.| 
[20071 [Devei-eux et al.|2009l >. 

Our best-fitting mass models include NFW haloes with a me- 
dian dark mass within the virial radius of = lO'^ '^^Mpj with 
rms scatter of 0.7 dex. With our adopted concentration-halo mass 
correlation, this corresponds to a concentration Cvir = 7.9 with rms 
scatter of 1.2. These parameters of the best-fitting dark haloes are 
defined in terms of their behaviour out to the virial radius, which is 
well beyond our kinematic constraints. They are therefore model 
dependent and should be treated with caution. The dark-to-total 
mass fraction within the galaxy is, however, a well-constrained 
quantity. We find that, on average, the haloes contribute around 15 
per cent by mass within and 49 per cent within /?25 • All but two 
galaxies are consistent with the maximal disk assumption defined in 
|Sackett|I997[ Models without dark matter are able to satisfactorily 
reproduce the observed kinematics in most cases, although there are 
problems at large radii in several galaxies. The improvement when 



a halo is added is statistically significant in all but four cases and the 
stellar mass-to-light ratios of mass models with dark haloes match 
the independent expectations of stellar population models better 

There is no systematic difference between the dark matter con- 
tent of the SO, Sa and Sb galaxies in the sample. This hints at a ho- 
mology between SOs and spirals. [Emsellem et al.|p007) and |Cap-| 
Ipellari, et al.| ( |2007 T| show that elliptical and SO galaxies exhibiting 
large-scale rotation (which they call fast-rotators) constitute a ho- 
mogeneous class in terms of their shape, stellar kinematics and pho- 
tometric properties. Assuming this direct link between SOs and fast- 
rotating ellipticals, our constraints on the dark matter content of Sa, 
Sb and SO galaxies can in turn be applied to fast-rotating ellipticals, 
implying a continuum of properties from Sb spirals to fast-rotating 
ellipticals. [Bertola et al.| l [T993^ made a particularly interesting pi- 
oneering suggestion when, with limited data, they suggested that 
spirals and ellipticals share a common scale, I.2/?e, within which 
luminous and dark matter contribute an equal amount to the total 
mass. For our larger sample, which contains no elliptical galaxies, 
we find that this radius is, on average, 4.1 /?e or 1.0i?25- 

The method we have used is generally applicable to axisym- 
metric (and spherically symmetric) galaxies and the code we have 
used is public. Although we have applied it to edge-on galaxies, it 
can equally be applied to any galaxy with a well-constrained incli- 
nation. We assumed that the galaxies are axisymmetric but argued 
that even in galaxies that are probably barred the method gives sen- 
sible results. Application of the method to simulated barred galaxies 
is needed to demonstrate this definitively, however Other limita- 
tions of our work include our inability to constrain anisotropy and 
the well-motivated but model-dependent way in which we eliminate 
one of the halo parameters. Kinematics above and below the major 
axis of the galaxy could allow us to lift these degeneracies without 
making assumptions, which would provide even more robust ob- 
servational tests of galaxy formation models and simulations. We 
have acquired long slit spectroscopy at multiple heights above the 
equatorial plane for a subset of the galaxies in sample, which will 
allow us to try this idea, but the ideal method is of course integral 
field spectroscopy from an instrument whose field of view and spa- 
tial resolution are optimized to reach the radii at which dark matter 
becomes important (e.g. |Verheijen et al.|2007^ . 



ACKNOWLEDGEMENTS 

We thank Aeree Chung for giving us access to the stellar kinematics 
of the sample and Giuseppe Aronica for giving us access to the Kn- 
band images. We also thank Liam Cook and Davor Krajnovic for 
participating in an early feasibility study for this project and Vic- 
tor Debattista, Richard Ellis, Susan Kassin and John Magorrian for 
valuable comments. We thank the anonymous referee for his/her 
helpful comments, which encouraged to consider the role of dark 
matter in more detail. MJW is supported by an STFC Postgrad- 
uate Studentship, MB by the STFC rolling grant 'Astrophysics at 
Oxford' (PP/EOOl 114/1) and MC by an STFC Advanced Fellow- 
ship (PP/D005574/1). This publication makes use of data products 
from the Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National 
Science Foundation. We also acknowledge the use of use of the 
HYPERLEDA database (|http : //leda . univ- lyonl . fr) and the 
NASA/IPAC Extragalactic Database (NED) which is operated by 
the Jet Propulsion Laboratory, California Institute of Technology, 



© 2009 RAS, MNRAS 000,[T||26l 



22 M. J. Williams et al. 



under contract with the National Aeronautics and Space Admin- 
istration. Numerical integration was performed using Craig Mark- 
wardt's QPINTlD numerical quadrature code. 



REFERENCES 

Abadi M. G., Navarro J. R, Fardal M., Babul A., Steinmetz M., 

2009, preprint (arXiv:0902.2477) 
Aguerri J. A. L., Debattista V. P., Corsini E. M., 2003, MNRAS, 

338, 465 

Arnouts S., et al., 2007, A&A, 476, 137 

Athanassoula E., Bureau M., 1999, ApJ, 522, 699 

AthanassoulaE., Misiriotis A., 2002, MNRAS, 330, 35 

Banerjee A., Jog C. J., 2008, ApJ, 685, 254 

Baugh C. M., 2006, Reports on Progress in Physics, 69, 3101 

Bedregal A. G., Aragon-Salamanca A., Merrifield M. R., 2006, 

MNRAS, 373, 1125 
Bell E. E, de Jong R. S., 2001, ApJ, 550, 212 
Bell E. E, Mcintosh D. H., Katz N., Weinberg M. D., 2003, ApJS, 

149, 289 
Bender R., 1990, A&A, 229, 441 
Bendinelli O., 1991, ApJ, 366, 599 

Benson A. J., Bower R. G., Erenk C. S., Lacey C. G., Baugh CM., 

Cole S., 2003, ApJ, 599, 38 
Bertola E, Pizzella A., Persic M., Salucci P, 1993, ApJ, 416, L45 
Binney J., Tremaine S., 2008, Galactic Dynamics: Second Edition. 

Princeton University Press 
Blanton M. R., Roweis S., 2007, AJ, 133, 734 
Blumenthal G. R., Eaber S. M., Elores R., Primack J. R., 1986, 

ApJ, 301, 27 

Blumenthal G. R., Eaber S. M., Primack J. R., Rees M. J., 1984, 

Nature, 311, 517 
Bolton A. S., Treu T, Koopmans L. V. E., Gavazzi R., Moustakas 

L. A., Buries S., Schlegel D. J., Wayth R., 2008, ApJ, 684, 248 
Borriello A., Salucci P, Danese L., 2003, MNRAS, 341, 1 109 
Bottinelli L., Gouguenheim L., Paturel G., Teerikorpi P., 1995, 

A&A, 296, 64 

Bullock J. S., Kolatt T. S., Sigad Y., Somerville R. S., Kravtsov 
A. v., Klypin A. A., Primack J. R., Dekel A., 2001, MNRAS, 
321,559 

Bureau M., Aronica G., Athanassoula E., Dettmar R.-J., Bosma 

A., Freeman K. C, 2006, MNRAS, 370, 753 
Bureau M., Athanassoula E., 1999, ApJ, 522, 686 
Bureau M., Athanassoula E., 2005, ApJ, 626, 159 
Bureau M., Ereeman K. C, 1999, AJ, 118, 126 
Burstein D., Davies R. L., Dressier A., Eaber S. M., Stone R. P. S., 

Lynden-Bell D., Terlevich R. J., Wegner G., 1987, ApJS, 64, 601 
Cappellari M., 2002, MNRAS, 333, 400 
Cappellari M., 2008, MNRAS, 390, 71 
Cappellari M., et al., 2006, MNRAS, 366, 1126 
Cappellari M., et al., 2007, MNRAS, 379, 418 
Cappellari M., McDermid R. M., 2005, Classical and Quantum 

Gravity, 22, 347 

Cappellari M., Neumayer N., Reunanen J., van der Werf P. P., de 
Zeeuw P T, Rix H.-W., 2009, MNRAS, 394, 660 

Chung A., Bureau M., 2004, AJ, 127, 3192 

Cole S., Aragon-Salamanca A., Erenk C. S., Navarro J. E, Zepf 
S. E., 1994, MNRAS, 271, 781 

Combes E, Debbasch E, Eriedli D., Pfenniger D., 1990, A&A, 
233, 82 

Combes E, Sanders R. H., 1981, A&A, 96, 164 



Croton D. J., Springel V., White S. D. M., De Lucia G., Erenk 
C. S., Gao L., Jenkins A., Kauffmann G., Navarro J. E, Yoshida 
N., 2006, MNRAS, 365, 11 
de Blok W. J. G., McGaugh S. S., 1997, MNRAS, 290, 533 
de Blok W. J. G., McGaugh S. S., Rubin V. C, 2001, AJ, 122, 
2396 

de Lorenzi E, Debattista V. P., Gerhard O., Sambhus N., 2007, 

MNRAS, 376,71 
de Souza R. E., Dos Anjos S., 1987, A&AS, 70, 465 
Debattista V. P, Sellwood J. A., 2000, ApJ, 543, 704 
Dekel A., Stoehr E, Mamon G. A., Cox T. J., Novak G. S., Pri- 
mack J. R., 2005, Nature, 437, 707 
Devereux N., Hriljac P., Willner S. P., Ashby M. L. N., Willmer 
C. N. A., 2009, in Jogee S., Hao L., Blanc G. Marinova I., eds, 
ASP Conf. Ser., Galaxy Evolution: Emerging Insights and Future 
Challenges, Astron. Soc. Pac, San Francisco, (arXiv:0902.0778) 
Devereux N. A., Becklin E. E., Scoville N., 1987, ApJ, 312, 529 
D'Onofrio M., Capaccioli M., Merluzzi P., Zaggia S., Boulesteix 
J., 1999, A&AS, 134, 437 
Eke V. R., Navarro J. E, Steinmetz M., 2001, ApJ, 554, 1 14 
Emsellem E., et al., 2007, MNRAS, 379, 401 
Emsellem E., Monnet G., Bacon R., 1994, A&A, 285, 723 
Englmaier P, Gerhard O., 1999, MNRAS, 304, 512 
Fioc M., Rocca-Volmerange B., 1997, A&A, 326, 950 
Fisher D., Illingworth G., Eranx M., 1994, AJ, 107, 160 
Freeman K. C, 1970, ApJ, 160, 811 

Gentile G., Salucci P, Klein U., Vergani D., Kalberla P, 2004, 

MNRAS, 351,903 
Gerhard O., Kronawitter A., Sagha R. P, Bender R., 2001, AJ, 

121, 1936 

Gerssen J., Kuijken K., Merrifield M. R., 2003, MNRAS, 345, 261 
Gnedin O. Y., Ki'avtsov A. V., Klypin A. A., Nagai D., 2004, ApJ, 
616, 16 

Jarrett T. H., Chester T., Cutri R., Schneider S., Skrutskie M., 

Huchra J. P, 2000, AJ, 1 19, 2498 
Jarrett T. H., Chester T, Cutri R., Schneider S. E., Huchra J. P, 

2003, AJ, 125, 525 
Jarvis B., 1987, AJ, 94, 30 
JarvisB. J., 1986, AJ, 91,65 
Jeans J. H., 1922, MNRAS, 82, 122 

Jensen J. B., Tonry J. L., Barris B. J., Thompson R. I., Liu M. C, 

Rieke M. J., Ajhar E. A., Blakeslee J. P, 2003, ApJ, 583, 712 
Jorgensen I., Eranx M., Kjaergaard P, 1996, MNRAS, 280, 167 
Karachentsev I. D., Karachentseva V. E., Parnovskij S. L., 1993, 

Astronomische Nachrichten, 314, 97 
Kassin S. A., de Jong R. S., Weiner B. J., 2006, ApJ, 643, 804 
Kauffmann G., White S. D. M., Guiderdoni B., 1993, MNRAS, 
264, 201 

Komatsu E., et al., 2009, ApJS, 180, 330 

Koopmans L. V. E., Treu T., Bolton A. S., Buries S., Moustakas 

L. A., 2006, ApJ, 649, 599 
Kranz T, Slyz A., Rix H.-W., 2003, ApJ, 586, 143 
Kuhlen M., Strigari L. E., Zentner A. R., Bullock J. S., Primack 

J. R., 2005, MNRAS, 357, 387 
Kuijken K., Merrifield M. R., 1995, ApJL, 443, L13 
Lutticke R., Dettmar R.-J., Pohlen M., 2000, A&AS, 145, 405 
Maccio A. V., Dutton A. A., van den Bosch F. C, 2008, MNRAS, 

391, 1940 
Maraston C, 2005, MNRAS, 362, 799 

Mayer L., Govemato E, Kaufmann T., 2008, Advanced Science 
Letters, 1, 7 



© 2009 RAS, MNRAS 000,[T||26l 



Constraints on stellar and dark matter in disk galaxies 23 



McGaugh S. S., de Blok W. J. G., Schombert J. M., Kuzio de 

Naray R., Kim J. H., 2007, ApJ, 659, 149 
Mei S., et al., 2007, ApJ, 655, 144 

Mendez-Abreu J., Corsini E. M., Debattista V. P., De Rijcke S., 

Aguerri J. A. L., Pizzella A., 2008, ApJL, 679, L73 
Merrifield M. R., Kuijken K., 1999, A&A, 345, L47 
Monnet G., Bacon R., Emsellem E., 1992, A&A, 253, 366 
Moster B. P., Somerville R. S., Maulbetsch C, van den Bosch 

F. C, Maccio A. V., Naab T, Oser L., 2009, preprint 
(arXiv:0903.4682) 

Mould J. R., Huchra J. P., Freedman W. L., Kennicutt Jr. R. C., 
Ferrarese L., Ford H. C., Gibson B. K., Graham J. A., Hughes 
S. M. G., lUingworth G. D., Kelson D. D., Macri L. M., Madore 
B. F, Sakai S., Sebo K. M., Silbermann N. A., Stetson P B., 
2000, ApJ, 529, 786 

Napolitano N. R., Capaccioli M., Romanowsky A. J., Douglas 
N. G., Merritield M. R., Kuijken K., Arnaboldi M., Gerhard O., 
Freeman K. C, 2005, MNRAS, 357, 691 

Navarro J. F, Frenk C. S., White S. D. M., 1997, ApJ, 490, 493 

Neto A. F, Gao L., Bett P, Cole S., Navarro J. F, Frenk C. S., 
White S. D. M., Springel V., Jenkins A., 2007, MNRAS, 381, 
1450 

Noordermeer E., Verheijen M. A. W., 2007, MNRAS, 381, 1463 

Palunas P, Williams T. B., 2000, AJ, 120, 2884 

Paturel G., Petit C., Prugniel P., Theureau G., Rousseau J., Brouty 

M., Dubois P, Cambresy L., 2003, A&A, 412, 45 
Persic M., Salucci P, Stel F, 1996, MNRAS, 281, 27 
Raha N., Sellwood J. A., James R. A., Kahn F. D., 1991, Nature, 

352,411 

Ratnam C., Salucci P, 2000, New Astronomy, 5, 427 

Roman N. G., de Vaucouleurs G., de Vaucouleurs A., Corwin Jr 

H. G., Buta R. J., Paturel G., Fouque P, 1991, Third Reference 

Catalogue of Bright Galaxies (RC3). Springer- Verlag 
Romanowsky A. J., Douglas N. G., Arnaboldi M., Kuijken K., 

Merrifield M. R., Napolitano N. R., Capaccioli M., Freeman 

K. C, 2003, Science, 301, 1696 
Rusin D., Kochanek C. S., Keeton C. R., 2003, ApJ, 595, 29 
Sackett P D., 1997, ApJ, 483, 103 

Schlegel D. J., Finkbeiner D. P, Davis M., 1998, ApJ, 500, 525 
Schwarzschild M., 1979, ApJ, 232, 236 
Scott N., et al., 2009, MNRAS, accepted. arXiv:0906.3321 
Shapiro K. L., Gerssen J., van der Marel R. P, 2003, AJ, 126, 2707 
Shaw M., 1993, A&A, 280, 33 

Shaw M., Wilkinson A., Carter D., 1993, A&A, 268, 511 
Shaw M. A., 1987, MNRAS, 229, 691 
Sheth K., et al., 2008, ApJ, 675, 1 141 

Silge J. D., Gebhardt K., Bergmann M., Richstone D., 2005, AJ, 
130, 406 

Skrutskie M. F, et al., 2006, AJ, 131, 1163 

Somerville R. S., Hopkins P. F, Cox T. J., Robertson B. E., Hern- 

quist L., 2008, MNRAS, 391, 481 
Springel V., White S. D. M., Jenkins A., Frenk C. S., Yoshida N., 

Gao L., NavaiTO J., Thacker R., Croton D., Helly J., Peacock 

J. A., Cole S., Thomas P., Couchman H., Evrard A., Colberg J., 

Pearce F, 2005, Nature, 435, 629 
Swaters R. A., 1999, PhD thesis, , Rijksuniversiteit Groningen, 

(1999) 

Swaters R. A., Madore B. F, van den Bosch F. C, Balcells M., 

2003, ApJ, 583, 732 
Thomas J., Jesseit R., Saglia R. P., Bender R., Burkert A., Corsini 

E. M., Gebhardt K., Magorrian J., Naab T., Thomas D., Wegner 

G. , 2009, MNRAS, 393, 641 



Thomas J., Saglia R. P., Bender R., Thomas D., Gebhardt K., 
Magorrian J., Corsini E. M., Wegner G., 2007, MNRAS, 382, 
657 

Tonry J. L., Dressier A., Blakeslee J. P., Ajhar E. A., Fletcher 
A. B., Luppino G. A., Metzger M. R., Moore C. B., 2001, ApJ, 
546, 681 

van Albada T. S., Sancisi R., 1986, Royal Society of London 
Philosophical Transactions Series A, 320, 447 

van den Bosch R. C. E., van de Ven G., Verolme E. K., Cappellari 
M., de Zeeuw P T, 2008, MNRAS, 385, 647 

Verheijen M. A. W., 1997, PhD thesis, PhD thesis, Univ. Gronin- 
gen, The Netherlands , (1997) 

Verheijen M. A. W., Bershady M. A., Swaters R. A., Andersen 
D. R., Westfall K. B., 2007, in de Jong R. S., ed. Island Universes 
- Structure and Evolution of Disk Galaxies, Astrophysics and 
Space Science Proceedings, Springer, Berlin, p. 95 

Wang L., Li C, Kauffmann G., De Lucia G., 2006, MNRAS, 371, 
537 

Wechsler R. H., Bullock J. S., Primack J. R., Kravtsov A. V, Dekel 

A., 2002, ApJ, 568, 52 
Weijmans A.-M., 2009, PhD thesis, Leiden University, (|http://| 

hdl.handle.net/1887/13970l 
Weijmans A.-M., Krajnovic D., van de Ven G., Oosterloo T. A., 

Morganti R., de Zeeuw P T, 2008, MNRAS, 383, 1343 
Weiner B. J., Sellwood J. A., WiUiams T. B., 2001, ApJ, 546, 931 
White S. D. M., Rees M. J., 1978, MNRAS, 183, 341 
Williams M. J., Bureau M., Cappellari M., 2009, in preparation 



APPENDIX A: MGE MODEL PARAMETERS 

For each galaxy, the parameters of the best-fitting MGE 
parametrizations of the projected light are presented in Table |Al1 



APPENDIX B: PHOTOMETRIC RECALIBRATION WITH 
THE 2MASS EXTENDED SOURCE CATALOG 

During the course of this work, we discovered that the photometric 
calibration zero points of |Bureau et al.| (j2006) were incorrect, so we 
recalibrated them using ^^-band images of the same objects taken 
from the Two Micron All Sky Survey (2MASS) Extended Source 
Catalog i |Skrutskie et al.|2006| >. 

We did this by measuring the light in matching elliptical aper- 
tures for each pair of images, and shifting the CASPIR image by a 
constant zero point offset such that its radial profile coincided with 
that of the 2MASS image. This offset gives a new, corrected A^^- 
band surface brightness zero point for each CASPIR image. 

Ensuring truly corresponding elliptical apertures enclosing 
identical locations on the sky is not trivial. The problem is that mea- 
surements of ellipticity and position angle based on two images of 
different depths will not necessarily yield identical results. Here the 
shallower 2MASS image does not reveal much of the faint discs, so 
important in determining the position angle and ellipticity. Corre- 
sponding apertures can trivially be ensured by using circular aper- 
tures. However, because circular apertures of increasing radii more 
quickly include noisy contributions from outside the galaxy (espe- 
cially in edge on systems), they do not use as much of either galaxy 
image as possible, decreasing the reliability of the recalibration. 

We therefore chose to impose the ellipticity measured for each 
CASPIR image on the corresponding 2MASS image and deter- 
mined the position angle independently for both images by using 
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Table Al. MGE parameters for the deconvoled ATs-band surface brightness. 
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Table Al. — continued 
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0.181 


+ 2.356 


1 719 


173 














IVrJP 71 9"^ 








IC5176 


















1 


4 812 


-0 289 


794 


1 


3.629 


-0.206 


0.746 
















+ 


A 1 77 


U.ZOo 


u. / 


+ 


3.480 


0.181 


0.746 


















3.884 


0.655 


0.634 




2.960 


0.466 


0.746 
















+ 


3.379 


0.905 


0.640 


+ 


2.979 


0.761 


0.640 


















4.335 


1.071 


0.114 


+ 


3.176 


1.126 


0.111 
















+ 


4.376 


1.107 


0.100 


+ 


3.346 


1.314 


0.103 
















+ 


2.215 


1.217 


1.000 


+ 


3.321 


1.333 


0.268 
















+ 


2.671 


1.224 


0.504 




3.803 


1.425 


0.181 
















+ 


3.069 


1.439 


0.100 


+ 


3.765 


1.453 


0.156 
















+ 


2.255 


1.764 


0.144 


+ 


2.242 


1.771 


0.165 

















Notes. Column (i) Sign of term in the Gaussian sum. (ii) Logarithm of the Gaussian amplitude, (iii) Logarithm of the Gaussian width, (iv) Axial ratio of the 
Gaussian. See equation j2j. 



an initial image truncated at the same approximate surface bright- [Bureau et al.| ( |2006l >, which were too faint, typically by around 0.8 

ness. In truncating both images, we temporarily remove information mag arcsec~^. 
from the deeper image so that it is missing as much of the disc as 
the shallower image. 



Before use, we correct the 2MASS images for the effects fore- 
ground Galactic extinction using the dust maps of [Schlegel et al.| 
(jT998l. This ensures that our recalibration incorporates this cor- 
rection. The calibration also involves an incidental and very slight 
colour transformation from Kg to A'n. Much like the filter used 
for 2MASS, the purpose of the /fn-band filter is to reject the ther- 
mal background admitted by a standard A'-band filter at its long 
wavelength end (Ski^utskie et al. 2006). We adopt a value for the 
absolute magnitude of the Sun at A'5-band of Mk^.q = 3.29 ( |Blan-| 
|ton & Roweis|200T] l. The result of this recalibration is presented in 
Table |Bl| The quantity A is the recalibration constant which should 
be added to the incorrectly calibrated A^n-band images presented by 
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Table Bl. Ks-hanA calibration corrections applied to the surface brightness 
of the images presented by |Bureau et al.|{2006} . 



Galaxy 



Galaxy 



(mag arcsec ) 



(mag arcsec ) 



B/PS huh 


^es 




Control 


sample 


NGC 128 


-0 


86 


NGC 1032 


-0.87 


ESO 151-G004 


-0 


80 


NGC 3957 


-0.73 


NGC 1381 


-0 


75 


NGC 4703 


-0.78 


NGC 1596 


-0 


84 


NGC 5084 


-0.97 


NGC 1886 


-0 


80 


NGC 7123 


-0.92 


NGC 2310 


-0 


96 


IC5176 


-0.82 


ESO311-G012 


-0 


83 






NGC 3203 


-0 


78 






NGC 3390 


-0 


83 






NGC 4469 


-0 


82 






NGC 4710 


-0 


85 






PGC 44931 


-1 


08 






ESO 443-G042 


-0 


75 






NGC 5746 


-0 


85 






IC 4767 


-0 


76 






NGC 6722 


-1 


29 






NGC 6771 


-0 


82 






ESO 185-G053 


-0 


70 






IC 4937 


-1 


21 






ESO 597-G036 


-0 


64 






IC 5096 


-0 


88 






ESO 240-GOll 


-0 


82 
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